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Abstract 

Several problems in mathematical physics relating to excitons in two dimen- 
sions are considered. First, a fascinating numerical result from a theoretical 
treatment of screened excitons stimulates a re-evaluation of the familiar two- 
dimensional hydrogen atom. Formulating the latter problem in momentum 
space leads to a new integral relation in terms of special functions, and fresh 
insights into the dynamical symmetry of the system are also obtained. A dis- 
cussion of an alternative potential to model screened excitons is given, and the 
variable phase method is used to compare bound-state energies and scattering 
phase shifts for this potential with those obtained using the two-dimensional 
analogue of the Yukawa potential. The second problem relates to excitons in 
a quantising magnetic field in the fractional quantum Hall regime. An exci- 
ton against the background of an incompressible quantum liquid is modelled 
as a few-particle neutral composite consisting of a positively-charged hole and 
several quasielectrons with fractional negative charge. A complete set of exci- 
ton basis functions is derived, and these functions are classified using a result 
from the theory of partitions. Some exact results are obtained for this complex 
few-particle problem. 

Introduction 

Dimensionality plays an important role in condensed matter physics, with many 
fascinating new phenomena emerging as the number of dimensions is reduced. With 
the introduction of improved semiconductor growth techniques over the past few 
decades, physical realisation of these systems has become possible. 

One of these systems is a quantum well, fabricated by sandwiching a layer of 
semiconductor between two layers of another material with a higher bandgap. The 
motion of electrons and holes is 'frozen' in the growth direction, and in modelling 
the system we may consider the electrons and holes as being essentially confined to 
a plane. As in bulk semiconductors, the Coulomb interaction between electrons and 
holes leads to bound states, known as excitons, which are extremely important for 
the optical properties of the quantum well. 
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The introduction of a strong magnetic field in the growth direction leads to 
further interesting phenomena in two-dimensional systems. At sufficiently low tem- 
peratures, many-electron correlations become important and new states of matter 
are formed. One of these states of matter is an incompressible quantum liquid, 
whose excitations carry fractional charge and obey fractional statistics. Spectro- 
scopic techniques can be used to probe the properties of such systems, and excitonic 
effects are again of the utmost importance. 

This Chapter will focus on a diverse range of problems in mathematical physics 
which have arisen as a result of our work on the theory of excitons in semiconductor 
nanostructures. The areas covered include the theory of special functions, dynamical 
symmetries in quantum mechanics, symmetric polynomials, the theory of partitions 
and counting of bound states in short-range potentials. 

In Section ^ we discuss a fascinating numerical result that arisies from a theo- 
retical treatment of screened excitons. This stimulates a re-evaluation of the famil- 
iar two-dimensional hydrogen atom. Formulating the latter problem in momentum 
space leads to a new integral relation in terms of special functions, and fresh insights 
into the dynamical symmetry of the problem are also obtained. A discussion of an 
alternative potential to model screened excitons is given, and bound-state energies 
and scattering phase shifts for this potential are compared with those obtained using 
the two-dimensional analogue of the Yukawa potential. 

In Section [21 we consider an exciton consisting of a valence hole and several 
fractionally-charged quasielectrons (anyons). A complete set of basis functions is 
obtained for this 'anyon exciton' and these functions are classified using a result 
from the theory of partitions. Expressions are derived for the overlap and interaction 
matrix elements of four- and six-particle systems, which describe an exciton against 
the background of incompressible quantum liquids with filling factors 1/3, 2/3 and 
2/5. Some exact results are obtained for the general (N + l)-particle case, including 
the binding energy for an exciton with zero in-plane momentum and zero internal 
angular momentum. 

1 Two-dimensional exciton: screened and unscreened 
1.1 Motivation 

Modelling of a screened Wannier-Mott exciton in two dimensions using the Stern- 
Howard potential X has revealed an intriguing numerical relation [21 Ej. In the 
search for an analytical derivation of this result, we address the simpler problem 
of an unscreened two-dimensional (2D) exciton, which has a well-known solution 
in real space. It has been claimed that the aforementioned numerical relation is 
not strictly exact, and an alternative approach to the problem has been proposed 
[3j. This has led to a new potential for modelling screened excitons, which is also 
considered here. 
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1.2 Screened 2D Wannier-Mott exciton 

Illumination of a semiconductor with near-bandgap light leads to the creation of 
electron-hole bound states known as excitons. The study of excitons in condensed 
matter physics has a long history 5 . Broadly speaking, excitons can be divided 
into two types: hydrogen atom-like Wannier-Mott excitons and more tightly-bound 
Frenkel excitons. In the present work we shall only be concerned with the first of 
these. 

Advances in semiconductor growth techniques have led to the creation of almost 
ideal 2D systems. In such systems, confinement leads to a higher exciton binding 
energy than in the bulk, which makes 2D excitons more experimentally accessible. 
A semiconductor quantum well is a quasi-two-dimensional system, in which pho- 
toexcited electrons and holes are essentially confined to a plane as their motion is 
'frozen' in the growth direction. Excitons are extremely important for the opti- 
cal properties of the quantum well. If enough electrons and holes are excited, the 
Coulomb attraction between particles may be weakened - this phenomenon is known 
as screening. The effect of screening on an exciton may be approximated by consid- 
ering the electron and hole as interacting via the so-called Stern-Howard potential 
PP, which has the following form in real space: 

V s {p) = -- jl - ^q s p[H. (q s p) - N (q s p)}} = --{1- q a p S Q , (q 8 p)} , (1) 

where q s is known as the 2D screening wavenumber, and Ho, Ao and So,o are the 
Struve, Neumann and Lommel functions, respectively. Note that all lengths and 
energies are measured in units of three-dimensional (3D) exciton Bohr radius and 
excitonic Rydberg, respectively. In these units, the radial equation for the relative 
electron-hole motion is given by 

d 2 R{p) IdR(p) / rl , m 2 \ , , __. . , n . 

7T T^+[Vs(p) + —)R(p) = ER(p). (2) 

dp z p dp \ p z J 

This is the radial part of the 2D Schrodinger equation, which is discussed in more 
detail in Sections 11.61 and 13.11 In Eq. (j2j), m is the azimuthal quantum number, 
which takes integer values. 

Numerical calculations Ej using Eq. Q have shown that with increasing 
screening, bound states in this potential disappear at integer values of the inverse 
screening wavenumber l/q s . This 'empirical' observation can be formulated as the 
following Sturm-Liouville problem: 

(A - F(x) + XV(x)F(x) = 0, (3) 
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where 



■t 1 d d 

A = -—x—, (4) 

x dx dx 



A = 2/q s , and we have made the substitution x = q s p. The zero right-hand side of 
Eq. (jHJ) reflects the fact that we restrict our interest to zero-energy bound states. 
The function V{x) is defined as 

V{x) = -- S 0>0 (x), (5) 
x 

and it has the following properties: 

V(x)xdx = l and V(x) = A ji - V(x)\ . (6) 



The key result of is that the solution F of Eq. (J2J) is square-integrable, i.e. 

poo 

/ \F(x)\ 2 xdx = 1, (7) 
Jo 

only if 

A = (2\m\ +v)(2\m\ +v + 1), (8) 

where v = 0, 1, 2, . . . is the number of non-zero nodes in the radial wavefunction. 

Notably, the potential defined in Eq. (^Q) is simplified considerably if we move to 
Fourier space: 

V s (q) = -7-^4- (9) 

W + Is) 

Eq. © is the 2D analogue of the familiar Yukawa potential (a standard textbook 
derivation for both the 2D and 3D cases can be found in [0], for example). 

The most compact way to reformulate the above Sturm-Liouville problem is in 
terms of a homogeneous integral equation: 



It follows that ^(q) is square-integrable in 2D Fourier space when A satisfies Eq. (jHJ). 
It should be emphasised that this result is of a numerical nature and has not been 
proved analytically. 
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In an attempt to find an analytical derivation of Eq. (jHJ) we address the simpler 
problem of an unscreened 2D exciton. This system is analogous to the 2D hydrogen 
atom, with the proton replaced by a valence hole, so we now focus our attention on 
such a system. 



1.3 Hydrogen atom and accidental degeneracy 

The in-plane motion of an electron and hole can be described by a two-dimensional 
Schrodinger equation for a single particle with a reduced mass. As in three dimen- 
sions, the 2D hydrogen atom exhibits 'accidental' degeneracy of the bound-state 
energy levels. 

Before proceeding it is necessary to discuss accidental degeneracy in more detail. 
The phenomenon of accidental degeneracy is well known from the theory of the 3D 
hydrogen atom. There are three quantum numbers {n,l,m}, although the energy 
only depends on the principal quantum number n. 

Degeneracy in m is of geometrical origin, and is due to the invariance of the 
system under spatial rotations. From Noether's theorem [7], this symmetry implies 
the existence of a conserved quantity, which is the angular momentum L. Indeed, 
quantum mechanics tells us that whenever there is a conserved quantity whose 
components cannot be measured simultaneously, then the energy eigenvalues will 
be degenerate [5|. 

The origin of the degeneracy in I, however, is not so obvious. In fact, this de- 
generacy does not have a geometrical origin, but arises from a hidden dynamical 
symmetry of the system. The conserved quantity associated with this accidental 
degeneracy is the Runge-Lenz vector A, the components of which cannot be mea- 
sured simultaneously. The existence of the Runge-Lenz vector is not so well known 
because when the Schrodinger equation is separated in spherical polar coordinates, 
the natural choice for a complete set of commuting operators is {H , L z , L 2 }, where 
H is the Hamiltonian. However, the presence of the Runge-Lenz vector is evident 
if the Schrodinger equation is separated in parabolic coordinates, as the natural 
choice for a complete set of commuting operators is then {H , L z , A z }. In general, 
accidental degeneracy appears when the Schrodinger equation can be separated in 
more than one coordinate system (for a detailed review see (Hj). 

In two dimensions there are two quantum numbers {n, m}. However, the energy 
eigenvalues only depend on n, so they are degenerate in the quantum number m. 
In this case the angular momentum only has a z-component (where z is in the 
direction perpendicular to the plane of motion), and so it can be known precisely. 
The degeneracy in m is therefore due to the existence of a 2D Runge-Lenz vector 
and a corresponding dynamical symmetry of the system. 
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Figure 1: Classical Runge-Lenz vector A. 



1.4 Historical background 

The Runge-Lenz vector has its origin in classical mechanics (see Fig. 0, and its 
discovery has been variously attributed to Laplace, Runge, Lenz and Hamilton (see 
jlUj and references therein). A result known as Bertrand's theorem states that the 
only central forces which result in closed orbits are the inverse square law (Kepler 
problem) and Hooke's law |1U1 111 j . The closure of orbits is due to the conservation 
of an additional quantity (apart from the energy and angular momentum); in the 
case of the Kepler problem this quantity is the Runge-Lenz vector 



where /i is the reduced mass, p is the momentum, L is the angular momentum, k is 
a positive constant and r is the radius vector of the orbit. 

A quantum mechanical analogue of the Runge-Lenz vector was first suggested 
by Pauli by analogy with classical mechanics. The main difference is that it is 
now an operator: 



where the form of the first term on the right-hand side ensures that A is Hermi- 
tian. Pauli also made the key conclusion that the existence of a Runge-Lenz vector 
constant in time implies that the hydrogen atom exhibits an additional degeneracy 



A = — (p x L) — k— 

// r 



(11) 




(12) 
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to that due to conservation of angular momentum. He then used A to derive the 
correct energy levels for the hydrogen atom. 

The most important paper concerning the symmetry of the hydrogen atom was 
published by Fock in 1935 13 . His first step was to consider the Schrodinger equa- 
tion in momentum space, which resulted in an integral equation. Considering only 
negative energy (bound-state) solutions, he then projected the 3D momentum space 
onto the surface of a four-dimensional (4D) hypersphere using a stereographic projec- 
tion. After a suitable transformation of the wavefunction, the resulting Schrodinger 
equation is seen to be invariant under rotations in 4D Euclidean space. The kernel 
is just the Green's function of Laplace's equation in four dimensions. Application 
of Green's theorem yields the correct energy eigenvalues, as well as eigenfunctions 
in terms of hyper spherical harmonics. Fock deduced that the dynamical symmetry 
of the hydrogen atom is described by the 4D rotation group SO (4), which contains 
the geometrical symmetry SO(3) as a subgroup. He related this hidden symmetry 
to the observed degeneracy of the energy eigenvalues. 

Shortly afterwards, Bargmann Jl] made the connection between Pauli's quan- 
tum mechanical Runge-Lenz vector and Fock's discovery of invariance under ro- 
tations in 4D space. The real-space Schrodinger equation was also separated in 
parabolic coordinates, which demonstrated more explicitly the conservation of the 
Runge-Lenz vector. 

A major step forward was taken by Alliluev JH]) who extended Fock's method 
of stereographic projection to the case of d dimensions (d ^ 2). He derived general 
expressions for the eigenfunctions and eigenvalues of the problem, and clarified some 
of Fock's arguments relating to the symmetry group for the d-dimensional case. 
Alliluev also considered the 2D harmonic oscillator, the symmetry group of which 
is isomorphic to the 3D rotation group SO(3). This provided an interesting parallel 
to the hydrogen atom case. 

The first explicit use of Fock's method for the 2D hydrogen atom was made 
by Shibuya and Wulfman [T2|. In this case, the eigenfunctions are just the usual 
spherical harmonics in three dimensions. The authors also included a graphical 
illustration of the degeneracy, in which two eigenfunctions which look very different 
on a plane are shown to be equivalent when projected onto a sphere. (It should be 
noted that the solution of the real-space Schrodinger equation in two dimensions 
was first provided in |17j . although the derivation was later repeated in |18j.) 

In the 1960s there was a revival of interest in accidental degeneracy as it was 
hoped that it could be applied to the theory of elementary particles. It attracted 
the attention of such famous names as Gell-Mann, Schwinger and Ne'eman. In a 
wide-ranging study |19| I20| , Bander and Itzykson provided a comprehensive review 
of work thus far. The first part gives a detailed explanation of the earlier work 
of Fock and Alliluev, the details of which are difficult to glean from the original 
works. Two systems of generalised spherical harmonics are then considered on the 
4D hypersphere: the standard hyperspherical harmonics, and an alternative in terms 



8 



D.G.W. ParBtt and M.E. Portnoi 



of Wigner D-functions. These are related to separation of variables is spherical polar 
and parabolic coordinates, respectively. For the 3D hydrogen atom, the Runge-Lenz 
vector in real space is shown to correspond to the infinitesimal generator of rotations 
on the 4D hypersphere. A detailed group-theoretical analysis of the hydrogen atom 
is also given. In the second part, the analysis is extended to scattering states, where 
the momentum space is now projected onto a higher-dimensional hyperboloid. 

The next few decades saw little progress, probably because the important ques- 
tions had already been addressed so thoroughly. However, improvements in semi- 
conductor growth techniques enabled the manufacture of effectively 2D structures, 
which soon led to a resurgence of interest in the 2D hydrogen atom. 

In a series of papers [22 [22l [22], the 2D Runge-Lenz vector was defined for 
the first time, and real-space solutions of the Schrddinger equation were applied to 
problems of atomic physics in two dimensions, including dipole matrix elements, the 
Stark effect, photon transition rates and fine/hyperfine structure. The relativistic 
case was also considered. 

Returning to higher dimensions, the use of hydrogenoid orbitals as basis sets for 
problems in quantum chemistry led to a re-examination of the separability of the 
(i-dimensional Schrddinger equation in different coordinate systems |24| . A novel ap- 
proach was taken by Nouri 25. , who considered the (i-dimensional Coulomb problem 
in quantum phase space. This led to an expression for the Wigner distribution func- 
tion, from which it is possible to calculate some physical and chemical properties for 
the hydrogen atom. The dynamical symmetry of the hydrogen atom has also been 
related to Kac-Moody algebras, rather than the usual Lie algebras |261 127j . This 
provides fertile ground for mathematical investigation and will no doubt lead to a 
new direction for research. 

Recent studies |28| I29| have shown that the classical Runge-Lenz vector has a 
deeper physical origin, being tied to the generator of Lorentz transformations for 
the relativistic two-body problem; the non-relativistic Kepler problem is then just 
a zero-order description of the more general case. In principle, it should be possible 
to undertake a similar analysis of the quantum mechanical problem based on the 
Dirac equation. 

Fock's method for the 3D hydrogen atom has been applied by Muljarov et al. |30[ 
131 j to the problem of an exciton with anisotropic mass and dielectric constant, 
providing an interesting alternative to a variational approach . The anisotropy 
is considered as a perturbation, and the perturbed wavefunction is expanded in 
terms of the hyper spherical harmonics. Brillouin- Wigner perturbation theory is 
then used to calculate the energy levels and eigenfunctions of the perturbed system. 

1.5 Fock's method 

We now present an outline of Fock's approach to the 3D hydrogen atom ^3], which 
led to important insights into the dynamical symmetry of the problem. We start 
from the Schrddinger equation in real space 
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-V 2 - - ) *(r) = EV(r), (13) 



r 



with distances and energies measured in terms of the 3D Bohr radius ao = fv 2 / fie 2 
and Rydberg Ry = n.e 4 /2h 2 , where n is the reduced mass of the system. 

Introducing the pair of Fourier transforms between real space and momentum 
space 



$(p) = J ^(r)e tpr dr, (14) 

* (r) = j^r I *( p ) e-ip ' r *' (15) 

and noting that for bound states the energy E = —p\ will be negative, we may write 
the Schrodinger equation as an integral equation 

The 3D momentum space is now projected onto the surface of a 4D unit hyper- 
sphere centred at the origin, and the in-plane momentum is scaled by po- A point 
on the hypersphere is specified by three angles, 6, 4> and a, where a is the angle 
between the position vector of the point and the fourth axis introduced in the 4D 
space. The transformation is given explicitly by 



u± = sin a sin V cos d> = k, (17) 

p 2 + po 

ii2 = sin a sin 9 sin <fi = ^°^ y _ ; (18) 



u$ = sin a cos V = », (19) 

p 2 + po 2 

P 2 ~ Po 2 , on , 

U4 = COS a = ~r, (20) 

p 2 + Po 2 



An element of surface area on the unit hypersphere is given by 

3 



(1p \ 
-r. o ) dp, (21) 
p 2 +Po 2 J 



and the distance between two points transforms as: 



10 



D.G.W. PaxGtt and M.E. Portnoi 



If the wavefunction on the sphere is expressed as 

, » 1 



xw 



2 i 2 \ 2 
P + P0 



/po v 2p 
then Eq. (|16|) reduces to the simple form: 

i r x(u') dn' 



*(p), 



(23) 



x( u ) 



2ir 2 p 



u — u 



'12 



(24) 



We now recognise that the kernel of the integral equation ()24|) is just the Green's 
function for the Laplacian operator in four dimensions: 



G 



u — u 



'12 ' 



which satisfies 



V 2 G = -4tt 2 <5(u - u') 



(25) 



(26) 



We also introduce a function tp which is a solution of Laplace's equation, i.e. it is 
a harmonic function. We now wish to apply Green's theorem for G and tp to the 
volume bounded by the unit hypersphere, letting n denote the unit vector from the 
origin to the north pole of the hypersphere: 



[tp(u')V 2 G - GV 2 ^(u')] dr 



dn dn 



dU. (27) 



However, this theorem can only be applied when both G and ip are harmonic inside 
the hypersphere and on its surface, whereas G has a singularity at u = u'. To 
overcome this problem we define our surface of integration to be the surface of the 
unit hypersphere with a small 'hemisphere' of radius e removed around the point u 
(see Fig.HJ): 



n e = {\u'\ 2 = 1, |u - u'\ 2 ^e}U {\u'\ 2 < 1, |u - u 

As e — > we regain the original unit hypersphere. 

Applying Green's theorem to this new surface we obtain 



,'|2 



(28) 



dn dn 



(29) 



because both G and ip are harmonic throughout this new region. Now, the integral 
on the right-hand side is made up of two parts: the first part will tend to the usual 
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Figure 2: Surface of integration in Green's formula. 



surface integral over the original unit hypersphere as e — » 0, whereas the second 
part over the small 'hemisphere' must be considered in more detail. 
For the first term, denoting the small hemispherical surface by u e : 



on 



1 3 



dS = 0(e) 



(30) 



where dS is an element of solid angle, and so this term vanishes as e — > 0. For the 
second term we have 



) as e -> 0. (31) 



Substituting back into Green's formula (|27|l and rearranging we obtain 



g Mh2 _ rt(u')— 

9r 9n 



tK2. 



(32) 



We now write our solution to Laplace's equation in the form 



V> = r A X , A = 0,1,2,. 



(33) 
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Substituting for ift in Eq. (|32j) and restricting u and v to the surface of the unit 
hypersphere, so that 

9 1 1 

7T 1 7V2 = ~ 1 7i2> 34 

on |u — u'| z |u — u'| z 

we obtain 

A + i r x (u')dn 
x(u) = - / l^W (35) 



Comparing this with Eq. (j2"I|) gives 

1 A + 1 

27T 2 J»o 27T 2 



(36) 



Substituting for an d associating A + 1 with the principal quantum number n, we 
find that 

E n = J, n= 1,2,3,... . (37) 



n- 



The corresponding eigenfunctions x( u ) are j us t the solutions of Laplace's equa- 
tion in four dimensions, i.e. the hyperspherical harmonics. These may be expressed 
alternatively as 



X (u) = (-2*ynJ 2n ^ 0! 1)! sin^C^^cosa)^^,^), (38) 



where Cjf are the Gegenbauer polynomials and Y™ are the conventional spher- 
ical harmonics. 

1.6 Two-dimensional hydrogen atom 

The relative in-plane motion of an electron and hole, with effective masses m e and 
nih, respectively, may be treated as that of a single particle with reduced mass 
/i = m e mh/(m e + mn) and energy E, moving in a Coulomb potential V(p). The 
wavefunction of the particle satisfies the stationary Schrodinger equation 



1 d ( d \ Id 



)2 



p dp \ dp ) p 2 dcf) 2 



*(p) = m( P ), (39) 



where (p, (ft) are plane polar coordinates. Note that excitonic Rydberg units are 
used throughout this Chapter, which leads to a potential of the form V(p) = —2/p- 
The eigenfunctions of Eq. (|39|) are derived in Section 13.11 It is well known that 
the bound state energy levels are of the form |17j : 
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E = n = 0,1,2,..., (40) 

(n+1/2) 2 V 7 

where n is the principal quantum number. Notably, Eq. H40J) does not contain 
explicitly the azimuthal quantum number m, which enters the radial equation (see 
Section|H3 Eq. (|241|) ). Each energy level is (2n + l)-fold degenerate, the so-called 
accidental degeneracy. The difference between the factor — l/(n + 1/2) 2 in the 2D 
case and — 1 /n 2 in the 3D case has been interpreted as arising from the topological 
peculiarity of 2D space [3I|: considering the origin as a singular point, paths that 
enclose the origin and those that do not cannot be continuously deformed into one 
another. 

It is convenient to introduce a vector operator corresponding to the z-projection 
of the angular momentum ~L Z = e z L z , where e z is a unit vector normal to the 
plane of motion of the electron and hole. We now introduce the 2D analogue of the 
quantum-mechanical Runge-Lenz vector as the dimensionless operator 

2 

A=(qxL z -L 2 xq) - -p, (41) 

P 

where q = — iV is the momentum operator. Note that A lies in the plane and has 
Cartesian components A x and A y . 

L z , A x and A y represent conserved quantities and therefore commute with the 
Hamiltonian: 

[H,L z ] = [H,A x ] = [H,A y } = 0. (42) 
They also satisfy the following commutation relations: 

[L z ,A x ]=iA y , (43) 
[L z ,A y ] = -iA x , (44) 
[A x ,A y ] = -AiL z H. (45) 

The existence of the non-commuting operators A x and A y , representing con- 
served physical quantities, implies that the Runge-Lenz vector is related to the 
accidental degeneracy of the energy levels in two dimensions (SJ. We now present a 
simple interpretation of the hidden symmetry underlying this degeneracy. 

For eigenfunctions of the Hamiltonian we can replace H by the energy E, and 
defining 

A' = — (46) 
2 V /Z £ 

we obtain the new commutation relations: 

[L z ,A' x ]=iA' y , (47) 
[L z ,A' y ] = -iA' x , (48) 
[A' x ,A' y ]=iL z . (49) 
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If we now construct a 3D vector operator 

J = A' + L 2 , (50) 

then the components of J satisfy the commutation rules of ordinary angular mo- 
mentum: 

[Jj, J fc ] = ie jM Ji, (51) 

where ejki is the Levi-Civita symbol. Indeed, the components of J are isomorphic 
to the generators of the 3D rotation group SO (3). 
Noting that A' • h z = L 2 • A' = 0, we have 

j 2 = (A' + L,) 2 = A /2 + L 2 , (52) 

where the operator J 2 has eigenvalues j(j + 1) and commutes with the Hamiltonian. 
The operator J 2 is known as a Casimir operator of the group consisting of the 
components of J. From Racah's theorem we find that for a Lie group of rank 
I, there exists a set of I Casimir operators which are functions of the generators and 
commute with every operator of the group. The present group is of rank 1 as none 
of the generators commutes with any other, which means that J 2 is unique. 

We now make use of a special expression relating A 2 and L 2 , the derivation of 
which is given in Section f3. 21 

A 2 = (4L 2 + 1) + 4. (53) 
Substituting in Eq. (]52(1 and again replacing H with E, we obtain 



J 



2 



1 

Te L 



£(4L* + l) + 4 +L 2 . (54) 



Because [ H, J 2 ] = 0, an eigenfunction of the Hamiltonian will also be an eigen- 
function of J 2 . Operating with both sides of Eq. (|54|) on an eigenfunction of the 
Hamiltonian, we obtain for the eigenvalues of J 2 : 

*7 + l) = -Q + i)- (55) 

Rearranging, and identifying j with the principal quantum number n, we obtain the 
correct expression for the energy eigenvalues: 

E = — ( * 2 . " = 0,1,2,.... (56) 
(n + 1/2) 



Note that the z-component of J is simply L z . Recalling that the eigenvalues of 
L z are denoted by m, there are (2j + 1) values of m for a given j. However, as j = n, 
we see that there are (2n + 1) values of m for a given energy, which corresponds to 
the observed (2n + l)-fold degeneracy. 
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Figure 3: Projection of 2D momentum space onto unit sphere. 



1.7 Fock's method in two dimensions 
1.7.1 Stereographic projection 

The method of Fock |13j . in which a 3D momentum space is projected onto the 
surface of a 4D hypersphere, may be applied to our 2D problem |361 137j . We begin 
by defining a pair of 2D Fourier transforms between real space and momentum space: 

$(q) = J *( P y* p dp, (57) 

We shall restrict our interest to bound states, and hence the energy E = —q$ will 
be negative. 

Substitution of Eq. (|58|) in Eq. i[39|l yields the following integral equation for 
*(q): 

= <»> 

The 2D momentum space is now projected onto the surface of a 3D unit sphere 
centred at the origin, and so it is natural to scale the in-plane momentum by qo 
(see Fig. El adapted from our paper [USD- Each point on a unit sphere is completely 
defined by two polar angles, 9 and (f>, and the Cartesian coordinates of a point on 
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the unit sphere are given by 



u x = sin 9 cos <f> = ^°^ x „ , (60) 
q 2 + qo 2 

(61) 



q 2 + qo 2 ' 



q 2 - qo 2 



u z = cos 6 = — j • (62) 



q + qo 

An element of surface area on the unit sphere is given by 



dil = sin clt) (U.> - ( o 2g0 \ dq, (63) 

\q 2 + qo 2 J 



and the distance between two points transforms as: 



(q 2 + qo 2 ) 1/2 (q' 2 + qo 2 ) 1/2 



If the wavefunction on the sphere is expressed as 



X(u) = -U(^M *(q). (05, 

v 7 * V 2 qo 



then Eq. (|59f) reduces to the simple form: 



i [ X (u r )dn> 

X(u) = g / -; Tp. (66) 

2irq J u-u' 



1.7.2 Expansion in spherical harmonics 

Any function on a sphere can be expressed in terms of spherical harmonics, so for 
x( u ) we have 

CO I 

*( U )=E E 4mlT(M), (67) 
1=0 m=-l 



where Yj m {6,(j)) are basically defined as in 
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where (cos 9) is an associated Legendre function as defined in [3""]. The constant 
c\m is an arbitrary 'phase factor'. As long as |c/ m | 2 = 1 we are free to choose q m , 
and for reasons which will become clear in Section [1.7. 31 we set 

Q m = H) H - (69) 
The kernel of the integral in Eq. (|66|) can also be expanded in this basis as |40| : 

1 °° - A 

1 1 A=0^=-A 

Substituting Eqs. (|6*7|) and (|70~|) into Eq. (|""oT) we have 



oo z 2 °° °° 'i '2 



EEwi^EE E E 

Z=0 m=— Z Zi=0 12=0 mi=— 1\ mi=— 12 



1 



2l 2 + 



-A hmi Y^ (9', 4>')Y^ (9, <j>)Y^\9\ 0') «BV. (71) 



We now make use of the orthogonality property of spherical harmonics to reduce 
Eq. ("""""J to the following: 

00/ ~ OO l\ -. 

>wr(M) = -E £ -j-^^ircM). ( 72 ) 

Z=0 m=— Z Zi=0mi=— Zi 

Multiplying both sides of Eq. ((72*)) by Y™'*{9,4>) and integrating over e?f2 gives 

A nr n' = 77; ; TT^nm') (^3) 

g (2n + 1) 

where we have again used the orthogonality relation for spherical harmonics. The 
final step is to rearrange for qo and identify the index n with the principal quantum 
number. This enables us to find an expression for the energy in excitonic Rydbergs: 

E = -q 2 Q = - 1 2 , 71 = 0,1,2,.... (74) 
(n + 1/2) 

This is seen to be identical to Eq. (|40jl. 

For a particular value of n, the general solution of Eq. (J""""|) can be expressed as 

n 

Xn (u)= ]T A nm Y™(9,<p). (75) 
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Each of the functions entering the sum in Eq. (|75jl satisfies Eq. (j66|) separately. So, 
for each value of n we have (2n+l) linearly-independent solutions, and this explains 
the observed (2n + l)-fold degeneracy. 

We are free to choose any linear combination of spherical harmonics for our 
eigenf unctions, but for convenience we simply choose 

Xnm(u) = A nm Y™{e,<j>). (76) 

If we also require our eigenfunctions to be normalised as follows: 

then Eq. ()76j) reduces to 

Xnm (u) =2KY™{9,(j ) ). (78) 

Applying the transformation in Eq. (|65[). we can obtain an explicit expression 
for the orthonormal eigenfunctions of Eq. ()59|): 

<Wq) = c wm j2 ^ n ~| m ||; (^^y 72 ^^^)^, (79) 
y (n+H)! W + QoJ 

where we have used the fact that qo = (n + 1/2) -1 , and 8 and <fi are defined by 

Eqs. dnnj-jnij). 

1.7.3 New integral relation 

To obtain the real-space eigenfunctions ^(p) we make an inverse Fourier transform: 

If i f'Z'K roc 

* (P) = Wfl $ ( q ) e " qPdq= (2^)2 J I H^e- iqpco ^'qdqd4>', (80) 

where cp' is the azimuthal angle between the vectors p and q. However, if we now 
substitute Eq. (|79[) into this expression we have to be careful with our notation. 
The angle labelled 4> in Eq. (|79*|) is actually related to 4>' via 

4> = <t>' + P , (81) 

where (ft p is the azimuthal angle of the vector p, which can be treated as constant 
for the purposes of our integration. 

Taking this into account, the substitution of Eq. ()79j) into Eq. ((50}) yields: 
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lPJ (27r)3/2\/ (n + |m|)! 



2-7T /■OO / r> \ 3/2 

" /n » (cos e)^^' cos ^ qdqdcf)'. (82) 



lo Jo W + Qo 
From Eq. (|62|) we obtain 

4 m l(cos^) = 4 m l/'^o!y (83) 

W + qo 2 J 

and we use the following form of Bessel's integral |4Uj : 

e i{m^- qp cos 4>') = 2ir(-i) m J m (qp), (84) 



where J m (qp) is a Bessel function of the first kind of order m. 
Substituting Eqs. (jE3J) and (jHU into Eq. (jH3) leads to 



, , = cnm(-i) m (n-\m\)\ im4>o 
{P) \(n+\m\)l 



oo 



2q 



o W + Qo 



3 /2 / 2 _ 2\ 

pM(£_tt Jm( ) d (85) 



We now make a change of variables, x = qop and y = q 2 /q$, so that Eq. (|85|) 
becomes 



*(p) = C r , 



I — ~\ \n+m(_-\m / Qo( n \m\)\ j m( j )o 

y 7r(n + \m\)\ 



X 







where we have used the fact that |4Uj : 

p\m\ ( V ~ l \ _ (_-i\n+mp\m\ I 1 ~ V \ (^\ 

n U+v n vi+i/y 

If we now equate the expression for ^(p) in Eq. Q86[) with that derived in Section 
13.11 we obtain the following: 

/■OO 
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The value of c nm chosen earlier in Eq. ensures that both sides of Eq. (|5H)) are 
numerically equal. For m ^ the relation simplifies to 

pm (i-y\ Jm(xyy) (-i)"(2x)- e - 2m 

n, m = 0, 1, 2, . . . ; m ^ n. (89) 

As far as we can ascertain, this integral relation between special functions has 
not previously been tabulated. For n, m = we recover the known integral relation 
1331: 

(l + yr /2 



1.8 Dynamical symmetry 
1.8.1 Infinitesimal generators 

Consider now a vector u from the origin to a point on the 3D unit sphere defined in 
Section ri.7.11 If this vector is rotated through an infinitesimal angle a in the (u x u z ) 
plane, we have a new vector 

u = u + 6u, (91) 
where the components of u are given in Eqs. (|60|) - (|62|) . and 

8u = a e y x u. (92) 

This rotation on the sphere corresponds to a change in the 2D momentum from q 
to q'. The Cartesian components of Eq. 1)91)1 are then found to be 

<=-^ = 4^+4^4 (93) 

q' 2 + q 2 q 2 +q 2 q 2 + ql 
, 1q qi y 2q q y 

U V q ,2 + q 2 q 2 +q 2> ^ 

n ' 2 = ?^ = ^^~ VT^' (95) 

where g 2 = q 2 + g 2 .. 

After some manipulation we can also find the components of Sq = q' — q: 

5q x = a q2 -f- 2q * , (96) 
2q K J 

6qy = _ a ME. (97) 

qo 

The corresponding change in ^(q) is given by 
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(q 2 + %) 3/2 \ 2 1o dq x q dq y J L J 
We can write this as 

<y*(q) = -^-a4*(q), (99) 
2<?0 

where the infinitesimal generator is given by 

A * = WTW 1 ( {e ~ "° ~ 2,i) ^ " 2&% 4) {e + * 2) " 2 ' (100> 

We now make use of the following operator expression in the momentum repre- 
sentation: 

p = e x x + e y y = iV q , (101) 

and the commutation relation 

[p, /(q)]=iVJ, (102) 

to derive a more compact expression for A x : 

A x = (q 2 -q%)x-2q x {q- p) -3iq x . (103) 

By considering an infinitesimal rotation in the (u y u z ) plane we can obtain a similar 
expression for A y . 

Ay = (q 2 - ql)y - 2q y {c i ■ p) - 3iq y . (104) 

These expressions operate on a particular energy eigenfunction with eigenvalue 
—q 2 . If we move the constant —q 2 to the right of the coordinate operators in 
Eqs. ()103|) and (|104[) . and replace it with the Hamiltonian in momentum space, 1ti: 

A x = q 2 x + xH - 2q x (q ■ p) - 3iq x , (105) 
Ay = q 2 y + yH- 2q y (q ■ p) - 3iq y , (106) 

then A x and A y can operate on any linear combination of eigenfunctions. 

1.8.2 Relation to Runge-Lenz vector 

Recall the definition of the 2D Runge-Lenz vector in real space: 

A=(qxL 2 -L 2 xq)- -p. (107) 

P 



22 



D.G.W. PaxGtt and M.E. Portnoi 



Using L z = p x q, and the following identity for the triple product of three vectors: 

ax (b x c) = (a-c)b- (a-b)c, (108) 
we can apply the commutation relation [p, q] = i to rewrite Eq. (|1U7|) in the form: 

A = q 2 p + p(V-^ - 2q(q • p) - 3iq. (109) 

If we now return to the expression for the real-space Hamiltonian in Eq. (J39|). it 
is apparent that we may substitute 

q 2 - - = H, (110) 
P 

in Eq. (tTU9l to yield 

A = q 2 p + pH - 2q(q • p) - 3iq. (Ill) 

Comparing this with Eqs. (|105|) and 1)1061) . it is evident that the two components 
of the Runge-Lenz vector in real space correspond to the generators of infinitesimal 
rotations in the (u x u z ) and (u y u z ) planes. 



1.9 Postscript 

Although projection into momentum space has proved useful for the unscreened case, 
it offers no further insights into the analytical derivation of Eq. Q . An alternative 
approach to this problem has recently been taken by Tanguy 4 . He claims, using 
the following argument, that the numerical result in Eq. (JSJ) for the critical screening 
lengths is not strictly exact. 

In a semi-classical approximation |41l I42j , the 2D formulation of energy quan- 
tisation for zero-energy bound states in the Stern-Howard potential Q may be 
expressed as 

^J- p {l -^q s p[Ho(qsP) ~ N (q s p)]} dp = (u + n. (112) 

Rearranging, the following expression for A = 2/q s is obtained in the semi-classical 
approximation : 




(113) 



where 

1 = I°° \F~ [H ° (n2) ■ n ° {u2)] du ' (n4) 



Exactly-solvable problems for two-dimensional excitons 



23 



and u = y/q s p- We would expect Eq. (|113j) to be accurate for y > 1. 

Tanguy asserts that for Eq. (JHJ) to be true we would expect 2Z = n, so as to 
obtain the correct leading order dependence of X sc on v 2 . He evaluates the integral 
(|114[) numerically and obtains 



21 = 3.14057414, (115) 

from which he concludes that A sc differs from A, and so Eq. (jSJ) is not exact. Using 
the fact that the potential Q behaves as — 2/p at small p and as —2/ (q 2 p 5 ) as 
p — > oo, Tanguy proposes another potential 

V(p) = - 2 2 , (H6) 

p0- + isp) 



which has the same asymptotic behaviour as potential lfl|). and has zero-energy 
bound states when A = 2/q s satisfies Eq. (JSJ). The corresponding zero-energy radial 
wavefunctions are of the form 



oc p' 



iFi -2\m\ - u, 2\m\ + v + 1; 2\m\ + 1; 



(2\m\+u){2\m\+v+l) 



(117) 



where 2-F1 is a hypergeometric function. 

The potential in Eq. (|116|) is clearly much simpler than the Stern-Howard po- 
tential and could certainly be used to model the screening of 2D excitons. To 
provide a comparison between the two potentials we use a 2D modification of the 
variable phase method (see Section f3.3l for the details) to calculate the bound-state 
energies and scattering phase shifts for Tanguy's potential l|llfij) . and compare these 
with similar calculations for the Stern-Howard potential. The results of our calcula- 
tions 03] are shown in Figs. and [3 Notably, for a wide range of values of q s there 
is at least ten per cent difference between numerically calculated binding energies for 
the two potentials, despite the fact that both potentials support the same number of 
bound states for the same value of q s . The difference in the scattering phase shifts 
is less pronounced. 



1.10 Conclusions 

We have shown that the accidental degeneracy in the energy eigenvalues of the 
2D Kepler problem may be explained by the existence of a planar analogue of the 
familiar 3D Runge-Lenz vector. By moving into momentum space and making a 
stereographic projection onto a 3D sphere, a new integral relation in terms of spe- 
cial functions has been obtained, which to our knowledge has not previously been 
tabulated. This approach has also yielded the momentum-space eigenfunctions for 
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1.0 




screening length 1/q 

Figure 4: Bound-state energies function of inverse screening wavenumber 

l/q s for the first three m = states in Stern-Howard (solid lines) and Tanguy 
(dotted lines) potentials. All values are normalised by the corresponding energies 
for the unscreened potential, E m ^ u {q s = 0) = — (m + v + 1/2) ~ 2 . 



K 3 




momentum k 

Figure 5: Scattering phase shifts function of the in-plane momentum k for 

l/q s = 5. Solid lines: Stern-Howard potential; dotted lines: Tanguy potential. 
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bound states, which will prove useful in future applications. We have also demon- 
strated explicitly that the components of the 2D Runge-Lenz vector in real space 
are intimately related to infinitesimal rotations in 3D momentum space. Although 
no further insight has been gained into the screened exciton problem, this is more 
than compensated for by the unexpected beauty of the unscreened case. 

With regard to the numerical expression (JBJ for critical screening lengths, it is 
still possible that it remains exactly true for both the Stern-Howard and Tanguy 
potentials. Tanguy's results are inconclusive: the semi-classical approximation re- 
mains an approximation, and it is possible that rounding errors in the numerical 
evaluation of the integral 1)114(1 lead to its deviation from tt. Indeed, it is interesting 
that the difference between the two potentials changes sign at q s ~ 1, which allows 
for the fact that Eq. (JHJ) may be simultaneously true for both potentials. 

The obvious extension of the 2D hydrogenic problem is to consider scattering 
states as well as bound states, in order to obtain a complete set of basis functions. 
However, rather than projection onto a sphere, this would require the 2D momentum 
space to be projected onto a two-sheeted hyperboloid. This has been done for 
the 3D hydrogen atom j20j, but the eigenfunctions are found to consist of two 
solutions, arising from the switching between sheets of the hyperboloid. Because of 
this complexity it is unlikely that the momentum-space eigenfunctions will produce 
any useful integral relations when compared with real-space solutions. 

2 Anyon exciton: exact solutions 

2.1 Motivation 

Intrinsic photoluminescence experiments have revealed interesting features in the 
optical spectra of incompressible quantum liquids (IQLs) at filling factors v = 1/3, 
2/3 and 2/5 [HJ 0H1 EE] • The anyon exciton model (AEM) was developed to provide 
a model of an exciton consisting of a valence hole and three quasielectrons (anyons) 
against the background of v = 1/3 and 2/3 IQLs (37], and has provided major 
insights into the role of electron-hole separation in determining their optical spectra. 
Recent developments in experimental techniques for studying such systems [481 149| 
have rekindled our interest in developing a more general form of the model which 
may be applied to higher fractions. To this end, we seek a generalised form of the 
AEM for a neutral exciton made up of a valence hole and an arbitrary number N of 
anyons. The particular case of N = 5 allows us to consider an exciton against the 
background of an IQL with filling factor v = 2/5. 

2.2 Quantum Hall effects 

The quantum Hall effect (QHE) was discovered in 1980 by von Klitzing et al. |5Uj . 
who studied an effectively 2D electron system at high magnetic field and low tem- 
perature. The Hall conductivity a xy exhibits plateaus at integer multiples of e 2 /h 
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(corresponding to an integer number of filled Landau levels), and is accompanied 
by a vanishing diagonal conductivity cr xx . This quantisation is independent of the 
shape of the specimen or presence of impurities and is accurate to one part in ten 
million, leading to its use as a standard of electrical resistance. 

Another significant breakthrough was made in 1982 by Tsui et al. 51_, who 
investigated higher-mobility semiconductor samples at increased magnetic fields and 
lower temperatures. Plateaus in a xy were observed at fractional values of e 2 /h, 
accompanied by minima in a xx . It soon became apparent that this phenomenon, 
labelled the fractional quantum Hall effect (FQHE), could not be explained in terms 
of a single-electron formulation; any theory would need to take into account the 
strongly-correlated behaviour of many electrons. (For a comprehensive review of 
the integer and fractional quantum Hall effects see |52| I53j.) 

In 1983, Laughlin showed [21] that the stability of the system at fractional fill- 
ing factor (number of filled Landau levels) is due to the condensation of the elec- 
tron system into an incompressible quantum liquid (IQL). By analogy with Jastrow 
wavefunctions in the theory of 4 He, Laughlin proposed a many-electron variational 
wavefunction for filling factor v = 1/m (m odd) to describe the ground state of the 
IQL: 

$ = nfe-^rll ex p(-Ni 2 /4) , (118) 

j<k I 

where Zi are the complex electron coordinates scaled in units of magnetic length, 
and m is an odd integer. He also noted the existence of an energy gap between the 
ground and excited states of the IQL, and predicted that excitations would be quasi- 
particles with fractional charge - quasielectrons and quasiholes. The first evidence 
for the existence of energy gaps was soon provided by magnetotransport experiments 
[53] and spectroscopic methods (57)1 EZ], and magnetotransport measurements also 
provided a signature of charge fractionalisation |58| I59| I6flj . 

An alternative to Laughlin's method was provided by Haldane |61j . who consid- 
ered a system based on a spherical geometry in which the electrons are confined to 
the surface of a sphere with a magnetic monopole at the centre. Numerical diagonal- 
isation of systems with up to eight electrons for v = 1/3 in a spherical geometry [H2] 
confirmed the accuracy of Laughlin's wavefunction in describing the ground state, as 
well as the presence of a finite energy gap between ground and excited states. This 
work was later extended to v = 1/5 with similar results [HH]- The use of Haldane's 
spherical geometry in numerical models of an exciton against the background of an 
IQL is considered in Section |2~H 

Further progress on the nature of the quasiparticle excitations was made by 
Halperin |64j . who developed pseudo- wavefunctions to describe low-lying excited 
states in terms of a small number of quasiparticles added to the IQL ground state. 
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These pseudo-wavefunctions are of the form 

N 

* = P(*k) II ^ - ^) Q II ex P {-N 2 /(4m)} , (H9) 

j<k 1=1 

where Z{ are the complex coordinates of the N quasiparticles, a is a statistical factor 
and P{zi) is a symmetric polynomial in the variables 2j. Halperin then used this 
method to construct the whole hierarchy of fractional states. From the quantisation 
rules that determine the allowed quasiparticle spacings he also determined that the 
quasiparticles are anyons obeying fractional statistics, i.e. the wavefunction changes 
by a complex phase factor when two particles are interchanged. (For more details 
on anyons and fractional statistics see jHSHOS]-) 

In addition to charged quasiparticles, the IQL also supports dispersive neutral 
excitations. The lowest branch of these neutral excitations is labelled a magnetoroton 
|67| I68j. and can be pictured as a quasielectron-quasihole pair. Strong evidence of 
the existence of these excitations has been provided by inelastic light scattering 
(Raman scattering) experiments [69] . 

Returning to Laughlin's wavefunction, a picture began to emerge of magnetic 
'vortices' generated by flux quanta, where the position of each vortex corresponds 
to a zero of the many-body wavefunction. The v = 1/3 FQHE, for example, could 
be interpreted as the attachment of three flux quanta to each electron, forming 
composite bosons (70]. Stability of the system at this filling factor then corresponds 
to a Bose condensate of these composite particles. 

The failure of this model to account for higher-order states with odd denomi- 
nator, as well as even-denominator states, led to the development by Jain of the 
composite fermion model |71|I72|, Vortices are again attached to electrons, but the 
fractional quantum Hall effect becomes an integer quantum Hall effect for the new 
composite particles. For example, in the v = 1/3 state each electron is consid- 
ered bound to two flux quanta, such that the composite particle feels an effective 
field of one flux quantum per particle. These composites can then be considered 
as fermions filling one Landau level. The composite fermion model is generally 
accepted as providing the best current description of higher-order FQHE states (al- 
though not everyone subscribes to this view [ZS] ) , as well as a method for generating 
many-particle wavefunctions and their excitations. 

One further development that should be mentioned is the recent direct obser- 
vation of fractionally-charged quasiparticles in 'shot noise' measurements |74| 175) . 
This noise arises from a current of individual quasielectrons with charge — e/3 tun- 
nelling through a narrow barrier (see [75] for more details). Whereas similar claims 
|77j have been questioned in the past, the results do look extremely promising, and 
a great deal of work is being carried out to verify these observations. 
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2.3 Photoluminescence experiments 

Although the earliest experimental work on the FQHE was based on magnetotrans- 
port experiments, these only allowed energies near the Fermi level to be probed. 
The introduction of spectroscopic techniques was a significant step as it enabled the 
whole spectrum of energies to be studied via contactless measurements. The ideal 
system for such optical experiments is a GaAs-AlGaAs heterostructure. Donors are 
placed in the layer of higher-bandgap material (AlGaAs) and the electrons associ- 
ated with these donors move towards the region of lower-bandgap material (GaAs). 
This process, known as modulation doping, creates a high-mobility two-dimensional 
electron gas (2DEG) in the layer of GaAs, which may be used as the basis for ex- 
perimental study of the IQL. However, the electric field that is used to create the 
2DEG also repels photoexcited holes, which necessitates adding a second barrier to 
restrict these holes. The physical system can be modelled by two parallel planes 
separated by a distance of h magnetic lengths. It has been found that when h = 
the system possesses a hidden symmetry and is insensitive to electron correlations 
|78l I79j (This is analogous to the familiar Kohn theorem, which states that the 
excitation spectrum of a translationally-invariant system is unaffected by electron- 
electron interactions [50].) As a result, the frequencies of all allowed transitions 
coincide exactly with those of an empty crystal. For this reason, charge asymmetric 
systems (with h ^ 0) are the most interesting for experimental observation. 

The phenomenon of photoluminescence also requires some explanation. Illu- 
mination of the semiconductor heterostructure produces non-equilibrium electrons 
and holes, which may then recombine radiatively via two different mechanisms. In 
extrinsic photoluminescence, electrons recombine with holes localised on acceptor 
dopants, introduced via <5-doping, whereas in intrinsic photoluminescence electrons 
recombine with 'free' holes. 

Experiments to study extrinsic emission provided the first optical measurement 
of FQHE energy gaps [22 Ell, an d a corresponding theory has been developed with 
good success |81| . Of most interest to the present work, however, are experiments on 
intrinsic PL. The first such experiments were by Heiman et al. |44| I45j . who studied 
high-mobility GaAs- AlGaAs multiple quantum wells at varying filling factor. In 
this real system the hidden symmetry is broken, thus allowing electron correlation 
effects to influence the optical spectra. Several intrinsic emission lines were observed 
in the photoluminescence spectrum for v < 1, and the most significant feature 
was a doublet near v = 2/3 which had a similar temperature dependence to the 
magnetoresistance associated with the IQL in transport experiments. The authors 
later studied a higher-mobility single quantum well, which led to the observation of 
a similar doublet at filling factors v = 1/3 and 2/5. The results for v = 2/3 were 
also reported by Turberfield et al. |46| . 

The only experimental work to study the effects of electron-hole separation on 
intrinsic photoluminescence was that by Plentz et al. |48| . Their apparatus allowed 
the 2DEG density to be changed by using applied gate voltages. This enabled the 
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effective electron-hole separation (in units of magnetic length) to be varied, while 
keeping the filling factor constant (for more experimental details see 82 ). The PL 
spectra obtained for v = 1/3 show a pronounced doublet structure. With increasing 
electron-hole separation, the energy splitting of the doublet decreases and the higher- 
energy peak dominates. Refinements to this technique have recently been reported 
by another group |83l I84j , and a similar method, which uses variations in the light 
intensity to change the 2DEG density, has also been developed [4"9"1 185]. 

These experimental methods have a great deal of potential for further investi- 
gation of spatially-separated systems. Although little work has been done in this 
area, it is hoped that the present study and subsequent work will rekindle interest 
in this direction. 

2.4 Theoretical studies 

Any attempt to model intrinsic PL must take into account low-energy excitonic 
bound states (due to electron-hole Coulomb attraction) against the background of an 
IQL. The first such work was by Apalkov et d. [13111111111131101, who numerically 
diagonalised systems comprising a small number of electrons together with a valence 
hole on the surface of a Haldane sphere. The electron-hole interaction was modified 
to take into account the layer separation h, which allowed both symmetric and 
asymmetric systems to be studied. Energy spectra were calculated for a range of 
values of h and exciton branches were classified in terms of the exciton internal 
angular momentum L. For h 1.5/, where / is the magnetic length, there was a 
single exciton branch with L = 0, whereas for h > 1.5/ new branches appeared, with 
higher values of L reaching down to become the ground state. The branches were 
labelled as either tight excitons or anyon excitons. For small h, the L = tight 
exciton forms the ground state, and as h is increased, anyon exciton branches form 
a succession of ground states and the tight excitons move to higher energies. 

Emission spectra were also calculated, and these also showed a strong dependence 
on h. For h < 0.5/ a single peak was observed, whereas for h > 0.5/ a second lower- 
energy peak appears. A tentative interpretation of these results says that the upper 
peak is from direct transitions for exciton momentum k = 0, whereas the lower 
peak corresponds to magnetoroton-assisted transitions with finite k. The double- 
peak structure shows a qualitative agreement with experiment |45l I82j . The results 
of Apalkov et al. were supported by a similar numerical study by Zang and Birman 
[9"T] , who also produced PL spectra for different h and obtained the familiar doublet 
structure. 

In parallel with finite-size calculations, an analytical model of an exciton against 
the background of an IQL was formulated by Portnoi and Rashba [4T | l92" l l9"3" l l94" | l95]. 
termed the anyon exciton model (AEM). This model is applicable for electron- hole 
separation h > 21, and thus complements the numerical approach. Indeed, remark- 
able agreement with finite-size results is obtained for h ~ 21. Initially, the simplest 
model of a hole with two semions (quasielectrons with charge — e/2) was considered, 



30 



D.G.W. PaxGtt and M.E. Portnoi 



but this was soon extended to a more-realistic four-particle exciton consisting of a 
hole and three quasielectrons with charge — e/3. The results of the AEM showed 
a multiple-branch energy spectrum, and a full set of exciton basis functions were 
obtained. More details of the formulation and results of the AEM can be found 
below. 

2.5 Review of anyon exciton model 
2.5.1 Formulation 

We present here a brief overview of the anyon exciton model (AEM) applied to the 
case of a four-particle exciton [17] . The model considers a neutral exciton consisting 
of a valence hole with charge +e and three quasielectrons (anyons) with charge —e/3 
and statistical factor a, subject to a strong perpendicular magnetic field H. The 
hole and anyons move in two different planes separated by a distance of h magnetic 
lengths. (Note that from now on we measure all distances in terms of the magnetic 
length I = (ch/eH) 1 / 2 unless otherwise specified.) 

This four-particle exciton will have a total of four degrees of freedom if all the 
particles are in the lowest Landau level. As the exciton is neutral we may assign 
it an in-plane momentum k, which absorbs two of these degrees of freedom. The 
remaining degrees of freedom result in internal quantum numbers and a multiple- 
branch energy spectrum. Indeed, the two internal quantum numbers are assigned 
the labels L and M, and enumerate exciton branches. Note that for k = 0, L is 
related to the projection of the angular momentum L z . 

To simplify the description of the exciton, the hole and anyons are considered 
as moving in the same plane, and the separation h can be introduced later as it 
does not appear in the exciton wavefunction. It is also convenient to move from the 
in-plane hole and anyon coordinates and to the new 2D coordinates 




(120) 



where jl = 12, 23 and 31, together with the complex coordinates Zji = Xji + iyji. 
In Eq. ()120() . R is the coordinate of the geometrical centre of charge of the exciton, 
p is the position of the centre of negative charge relative to the hole, and rji are 
difference coordinates between pairs of anyons. The following constraint should also 
be noted: 

ria + r 23 + r 3 i = 0, (121) 

which greatly simplifies the subsequent formulation. 

The motion of an anyon exciton is described by a Halperin pseudo-wavefunction 
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*(R, p, {%,}) = -jL= exp |ik • R + £z • [p x R] - i(p- d) 2 | 

x Pl(. . . , z- jU . . .) exp {-|z j7 | 2 /36} , (122) 

jl 

where A is the area of the sample and Pl is a homogeneous polynomial of degree 
L in the difference coordinates. Note that the statistical factor a = for bosons, 
a = 1 for fermions, a = —1/3 for quasielectrons in a v = 1/3 IQL and a = +1/3 
for quasielectrons in a i/ = 2/3 IQL. 

The properties of the polynomial Pl are of the utmost importance. By definition 
it is symmetric in the anyon coordinates Zj, but the permutation rules with respect 
to the difference coordinates Zji are not obvious. Indeed, it turns out that if L is 
even, Pl is symmetric in Zji, whereas if L is odd, the polynomial is antisymmetric in 
Zji. For L = 3, the simplest case, this polynomial is a Vandermonde determinant: 



W (Z12,Z 2 3,Z 31 ) 
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Zl2 


Z23 


231 


7 2 


7 2 


7 2 


z 12 


z 23 


z 31 



(Zu ~ Z 2 3)(Z23 ~ Z3l)(^31 ~ ^n)- (123) 



The linearly-independent L-even polynomials can be chosen as 

p r i- - - \ _ -L-4M -2M -2M , -L-4M -2M -2M , -L-4M-2M-2M 94 x 
-TL,M (^12, ^23, 2 3 i J — Z 12 Z 23 Z 31 + Z 23 Z 31 Z 12 + Z 31 Z 12 Z 23 , (124J 

where the quantum number M enumerates different polynomials with the same 
degree L. The total number of L-even polynomials with a given L is therefore 
[L/6] + 1, where the brackets [ ] denote the integer part. The L-odd polynomials 
are constructed as follows: 



P 



3,0 



W, P L . 



M 



WP L - 



•3,Af i 



(125) 



and there are \{L — 3)/6] + 1 of these polynomials for a given L. 

Given an expression for the exciton basis functions *$>i in terms of the quantum 
numbers k, L and M, any exciton wavefunction may be expressed in terms of this 
basis: 



E 



(126) 



However, the basis functions are not orthogonal in M, so the Schrodinger equation 
is written as follows: 

H X = eB X , (127) 

where B is the overlap matrix of scalar products and e is an energy eigenvalue. The 
evaluation of the matrix elements H and B is rather complicated, and the reader 
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is referred to [I7j for further details. These matrix elements allow the calculation 
of the energy spectra and electron density for the exciton, and these results are 
outlined below. 

2.5.2 Results 

For energy eigenvalue calculations at k = 0, a comparison of the energy for statistical 
factor a = —1/3 and a = shows good agreement for h > 2 (see Fig. 1 in 4Tj ) , and 
hence a boson approximation is made to simplify further calculations. For k ^ 0, 
the energy spectra must be calculated numerically, and it is found that in the region 
h > 2 the ground state of the system has non-zero angular momentum L = 3m, 
where m ^ 2 is an integer. Note also that the wavefunctions of these ground 
states contain the polynomials Pqm,m and WPqm,m- These satisfy the hard-core 
constraint, in that they go to zero if any two of the anyon coordinates are the same. 
This is very encouraging as hard-core functions are the accepted standard for a 
proper description of anyons (HHl E7J • 

Energy spectra for h = 2 and 3 are shown in Fig. 37J- For h = 2, the 
minimum of the spectrum occurs for non-zero k, which means that the exciton in 
the lowest energy state possesses a non-zero dipole moment. The negative dispersion 
in Fig. Efa) arises because of the mutual repulsion of L = 2 and L = 3 branches. 
For higher values of h the ground state is always at k = 0, as in Fig. Efb) , and the 
lowest-branch dispersion is always positive. 

Calculations of the electron density around the hole also lead to remarkable re- 
sults (see Fig. [7J 47 j). For k = there is a crater-like dip in the density around 
the origin; this is an obvious signature of charge fractionalisation as a normal mag- 
netoexciton always has a maximum at this point. The distribution for non-zero k 
clearly shows two anyons close to the hole and a single anyon at a slight distance from 
it, providing further clear evidence of charge fractionalisation in anyon excitons. 

2.6 Exciton basis functions 
2.6.1 Derivation 

The concept of exciton basis functions is now generalised to an exciton consisting 
of a valence hole with charge +e and N anyons with charge —e/N and statistical 
factor a. The hole and anyons reside in two different layers, separated by a distance 
of h magnetic lengths, and are subject to a magnetic field H = Hz perpendicular to 
their planes of confinement. Unless otherwise stated, we shall assume that the hole 
and the quasielectrons are in their corresponding lowest Landau levels. To simplify 
the description of the exciton, the hole and anyons may be considered as moving in 
the same plane with coordinates and ij, respectively, and the the layer separation 
h can be introduced later when considering the anyon-hole interaction. 

An exciton consisting of a hole and N anyons, all in the lowest Landau level, 
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Figure 6: Anyon exciton dispersion e(k) for two values of electron-hole separation 
h. Numbers show L values; h is in units of magnetic length I and e is in units of 
the Coulomb energy e c = e 2 /el. 




Figure 7: Electron density distribution in an anyon exciton for different values 
of the exciton in-plane momentum k. The distance h between the hole and the 
incompressible electron liquid is equal to two magnetic lengths. The hole is at the 
origin; the x-axis is chosen along the exciton dipole moment. 
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will have a total of N + 1 degrees of freedom. Assigning the exciton an in-plane 
momentum k absorbs two of these degrees of freedom, so for N ^ 2 the exciton 
will have N — 1 internal degrees of freedom. At this stage we consider the particles 
as non-interacting and introduce interactions later as necessary. We can therefore 
write the (N + l)-particle Hamiltonian as 

where qh = +e and q a = —e/N are the hole and anyon charges, respectively. 
Choosing the symmetric gauge 

A = \ [H x r] , (129) 

and scaling all distances with the magnetic length / = (chj 'eH) 1 / 2 , we obtain 

#o = ( iv h - [z x r fc ]) J-Qvj + ipxrj]) , (130) 



2m h 



where, as usual, e, h and c are assumed equal to unity. 

We now introduce the following new coordinates (see Fig. |SJ) : 

1 / 1 N \ 1 N 1 N 

R= 2 r/l+ iV^ rj ' P = rh ~N^ rj > ^ =r ' _ ]vE r '' (131) 
\ i=i / j=i i=i 

together with the complex coordinates Q = £ x j + i£ y j. Note also the following 
constraint on these coordinates: 

N N 

E ^ = E = 0. (132) 



The fact that the variables £„■ are not independent means that derivatives with 
respect to these variables are not defined. However, it is possible to introduce the 
derivatives which treat the variables £ 3 - as if they were independent |99| . 

Let S be the 2A r -dimensional space of physical coordinates rj, and Sg be the 
2A^-dimensional space of the independent coordinates £„-. Due to the constraint 
above, physical motion only takes place within a submanifold Mg C 5^, so ^ maps 
the physical space S to the submanifold Mg, i.e. £ : S* — > Mg. An arbitrary complex 
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Figure 8: Two-dimensional coordinate system for hole and N anyons (C indicates 
the centre of negative charge). The hole is considered as being in the same plane as 
the anyons, with layer separation introduced as required. 



function / can always be expressed explicitly as a function of the coordinates £j , and 
so we can write that / maps Mg to C, i.e. / : — ► C. If we now keep the expression 
for / fixed, we can extend its definition over the whole space by removing the 
constraint, and it is now correct to talk about its derivatives V$./(. . . , £j, . . .). 

We may now consider a composition of the two maps S — ► C —* C (see 
Fig. IHJ). Wavefunctions of the kind used to describe motion with respect to the 
centre of charge of the anyon subsystem are compositions of / and £, such that 

(/ o £)(n, . . . , v N ) = /(... ,^(n, . . . , TN ), . . .). (133) 

We do not distinguish between / and / o £ in what follows, but rather write both 
as /. 

Using the chain rule, we may now rewrite the derivatives of / with respect to ta 

as 

N 

i=i 

1 - 

V € ./i.... -Y. V £ /(....£,....;. (134) 

i=i 

We may now express the old coordinates in terms of the new ones as follows: 

r fc = R+|, Tj =R-| + ^., (135) 
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Figure 9: Wavefunction of relative motion expressed as a composition of two map- 
pings S-^M^cS^C. 



and the derivatives with respect to the old coordinates are therefore 



v ft = -v R + v p , 



11 1 - 

V,- = V R Vn + Vt . 

i=i 



(136) 
(137) 



It is now possible to write the Hamiltonian (J13U|) in terms of the new variables. 
Using Eqs. (fTM |) -(fT37 )) we obtain 



Hn — H exc + He 



(138) 



where 



Hexc = { tttVr + - V p - - 

zrrih \2i i 2 L 



z x ( R 



1 f 1 „ !^ 1 



z x R + 



, (139) 



and 
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(140) 



Note that Eq. (|14U|) was obtained by applying constraint (|13'2|) . 

The first part H exc is similar to that for a standard diamagnetic exciton jlUUI 
I1U1) . with the electron mass replaced with N anyon masses. The eigenfunctions of 
this operator can be written straightforwardly as 

1W(R, p) = exp {ik • R + iz ■ [R x p] /2} (141) 

where <&(p) satisfies the equation 



i Vo _ [A ><i^Y + j_ 1 1 Vp + ^<(p^| 2 j tt = ^ 

(142) 



2m h \ i p 2 J 2iVm a 



Here, d = k x z is the exciton dipole moment. Note that Eq. (|141|) does not contain 
the mass-dependent phase factor that appears in the wavefunction in jlOlj . This is 
due to our choice of the centre of charge of the exciton for our coordinate R, rather 
than the centre of mass chosen in |100[ll0f] . Thus, the quantum number k entering 
Eq. 1)141)1 represents the momentum of the geometrical centre of the exciton. 
The solutions of Eq. 1)142)) are 

$ nm = | p |HlM ( (p _ d) 2 /2 } exp { imxj) _ (p _ d)2/4 } (143) 

where (p is the azimuthal angle of the vector (p — d) and L™(z) are the associated 
Laguerre polynomials. The corresponding eigenvalues are 

„ 1 / \m\ — m + l\ 1 / |m|+m + l\ , 
E exc = — n + J— ! + n + i— ! . 144) 

Note that in Eqs. (j!43j) and ()144)l we no longer restricted ourselves to the lowest 
Landau level for the hole and quasielectrons. For the ground-state (n,m = 0), 
Eq. 1)141)1 reduces to 

^ exc (R, p) = exp {ik • R + iz • [R x p] /2 - (p - d) 2 /4} , (145) 

which depends on the quantum number k alone. 

We now move to consider the term in Eq. (|14()j) . which may be divided into 
two parts: 



H 



1£ 



4) £(7 v «. + d?[* *«'])'■ (146> 

/ 1=1 x ' 
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and 



+ ^^l 2 +^,(gV Sl )-i(gV 5 ,) 2 }. (147) 

The eigenfunctions of the ground state of .ffij can then be written in terms of 
complex coordinates as 

*e = II(0) /3 ' ex p{-y 2 /4N} ) (148) 

i 

with the corresponding eigenvalues 

Ec = — ( 1 - — V (149) 

If we now apply constraint ()132|) , we find that the function is also an eigenfunction 
of H2£, with a corresponding eigenvalue of zero, i.e. 



. 2 

AT 



H^oc [Y J ^ \ ^C = 0. (150) 

A calculation of the total energy of the ground state gives 

E = E exc + E^ = -^— + -^— = l(u; h + Nuj a ), (151) 
znih 2m a 2 

which is indeed the ground-state energy of the original (N + l)-particle Hamiltonian 
in Eq. OB . 

The most general form for the ground-state eigenfunction of Hq is 

(R, p, {Ci}) = exp {ik ■ R + iz • [R x p] /2 - (p - d) 2 /4} 

x F(..., £,...) J] exp{-|C P | 2 /4iV}, (152) 
p 

where F is a function of the complex conjugates Q alone. We choose F as follows 
to satisfy the interchange rules for anyons: 
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F(. . . , Q, ...) = P L (..., . . .) H(Q - (if, (153) 

3<l 

where Pl is a symmetric polynomial of degree L in the variables Q. Note that 
for k = the problem has rotational symmetry about the z-axis and the de- 
gree of the symmetric polynomial L is related to the exciton angular momentum 
[L z = —L — N(N — l)a/2]. 

We are now in a position to express the anyon exciton basis functions in the final 
form 

(R, p, {&}) = exp {ik ■ R + iz ■ [R x p] /2 - (p - df/A] 

x P L (...,C i ,...)II(0-^) Q II ex P{-|^| 2 / 4Ar }- (154) 

j<l v 

2.6.2 Symmetric polynomials 

We now consider the precise structure of the symmetric polynomials which ap- 
pear in Eq. (|154|) . To determine the symmetric basis polynomials of order L, we 
apply the fundamental theorem of symmetric polynomials |l()2j . which states that 
any symmetric polynomial in N variables, P(x\, . . . , xjv), can be uniquely expressed 
in terms of the elementary symmetric polynomials 

cri = ^Xi, CT 2 = ^XiXj, a 3 = y~] XjXjX k , . . . , a N = X\x% ■ ■ ■ x N . (155) 

i i<j i<j<k 

It is apparent that all linearly-independent symmetric basis polynomials of a partic- 
ular degree may be enumerated by considering the possible products of the elemen- 
tary symmetric polynomials ()155|) . as the total degree L is the sum of the degrees 
of the constituent polynomials. For example, the possible symmetric polynomials 
of degree four are constructed from af, erf 02 > fios, u\ and 04. 

There are two further problems. The first is to determine the number of possible 
products of elementary symmetric polynomials for a total degree L. The second is 
to determine the structure of these products. 

To calculate the number of linearly-independent symmetric basis polynomials 
of degree L which may be constructed from the elementary symmetric polynomials 
<7l, . . . , (Tjv, we use a result from the theory of partitions jlU3j . 

The number of ways of partitioning a number L into parts of size 1,2, ... ,P is 
given by the coefficient of x in the expansion of 
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where each term in the product is known as a generating function. The num- 
ber of ways of partitioning increases rapidly with L and does not follow any pat- 
tern. Furthermore, the different products of elementary symmetric polynomials for 
a particular value of L must be determined by hand. An approximate formula for 
calculating the number of ways of partitioning a number does exist, the so-called 
Hardy- Ramanuj an formula jl04j . but it is not applicable to the present case as we 
have no polynomial of order one due to constraint (|132j) . This constraint signifi- 
cantly reduces the number of possible symmetric polynomials by removing the first 
factor in the product (j!56|) . For example, for N = 3 the number of polynomials of 
degree L is equal to the integer part of L/6 + 1 for even L and the integer part of 
(L — 3)/6 + 1 for odd L, which corresponds to the result obtained in [4*7] . 

2.6.3 Discussion 

The key difference between the current general formulation and that outlined in 
|47j for the four-particle case (N = 3) is the replacement of anyon difference co- 
ordinates Zji by the new coordinates Q. The principal advantage of the difference 
coordinates was that they simplified the calculation of inter-anyon repulsion matrix 
elements, as well as the form of the statistical factor in the exciton wavefunction. 
However, they had the disadvantage that the classification of symmetric polyno- 
mials was more difficult, as it was necessary to introduce a Vandermonde deter- 
minant for odd-L polynomials. For example, in the simplest case L = 3, even 
though the number of anyon coordinates is the same in both formulations, we have 
o" 3 = {z\2 — Z23){z23 — ^3i)(^3i — £12) in terms of difference coordinates, whereas we 
have a simple product o~ 3 = C1C2C3 i n terms of the new coordinates. For a number of 
anyons greater than three this disadvantage is crucial, as the classification of poly- 
nomials in terms of Zji becomes too cumbersome and the number of constraints on 
these coordinates is greater than one. 

2.7 Reformulation of four-particle anyon exciton 
2.7.1 Formulation 

We now use the above results to formulate the problem of a four-particle exciton 
in terms of the new coordinates Q. We consider an exciton consisting of a hole 
and three anyons with charge — e/3, and make a boson approximation so that the 
statistical factor a = 0. As mentioned in Section 12.5.21 it is evident from Fig. 1 
of |47j that for large values of h (required for the AEM to be valid) the statistical 
factor a becomes unimportant, and the results for a = were very similar to those 
for a = ±1/3. We are therefore justified in approximating the anyons as bosons. 
The basis functions for a four-particle anyon exciton in a boson approximation 
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Table 1: Possible ways of constructing a symmetric polynomial Pl from the elemen- 
tary symmetric polynomials 0"2 and cr 3 . 



Order, L 


No. of polynomials 


Structure 
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1 
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1 
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0"3 
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,1 
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O2O3 
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1 2 
*2. CT 3 
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<?2 0-3 
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4 2 
of, 0- 2 cr| 
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0"2°"3 5 0"3 


10 
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1 2 2 
»2 > ^2^3 


11 
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4 3 
0"2°"3, O"20"3 


12 


3 


0-2 > 02*3 > *3 



are, from Eq. I|154[) : 

*(R,p,Ci,C2,C3) =exp{ik-R + iz- [Rx p] /2 - (p - d) 2 /4} 

3 

xPL,A/(Ci,C2,C 3 )n ex P{-|^| 2 /12}, (157) 

p=l 

where the index M enumerates different linearly-independent symmetric polynomi- 
als of degree L. 

As a result of constraint ()132j) . o\ = Cl + C2 + (3 = 0? anci so we need to 
consider only the elementary symmetric polynomials 02 and 03 when constructing 
the symmetric polynomial Pl,m ■ From Eq. Q156|) . the number of possible ways of 
constructing a polynomial of degree L is the coefficient of x in the expansion 

1 1 

z = 1 + x 2 + x 3 + x 4 + x 5 + 2x 6 + x 7 + 2x 8 + . . . . (158) 

1 — x A 1 — X 6 

Tabled shows the possible ways of constructing the first twelve polynomials Pl. 
As in Section r2.5.11 a general exciton wavefunction for given k may be expanded 
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in terms of the complete set of basis functions Q157|) as follows: 

* = ( l59 ) 

i 

Basis functions with different values of k are orthogonal. We shall later show that 
functions with different L are also orthogonal. However, basis functions with the 
same value of L but different M are not necessarily orthogonal, so that their scalar 
products (L, M\L, M') will be non-zero. We therefore write the Schrddinger equation 
in matrix form as 

H X = eB X , (160) 
where B is the block diagonal matrix of scalar products (the overlap matrix) and 
e is an energy eigenvalue. The size of each block in B depends on the number of 
different wavefunctions with given L. In the case of a four-particle exciton (see 
Table the block corresponding to L = 12 will be of size 3x3. Both the overlap 
matrix B and the Hamiltonian matrix H are diagonal in k. Note that for k = the 
matrix H takes the same block diagonal form as B (as will be shown later), and as 
a result the problem becomes exactly solvable. 

We shall now proceed to evaluate the matrix elements in Eq. ()16U|) . As H and 
B are diagonal in k, in what follows we consider only the matrix elements diagonal 
in k. Since all terms in the polynomial Pl M are °f the form of a product of the 
coordinates (j, we use the following functions in monomials in the matrix element 
calculations: 

1- Wik (R, p, {£}) = C exp {ik • R + iz • [R x p] /2 - (p - d) 2 /4} 

x CrC 2 n2 C 3 n3 fl ex P HC//12} , (161) 
i=i 

where {n} denotes the set of quantum numbers m, ni and n^, and the constant 
C will be defined in the next Section. The basis functions in terms of symmetric 
polynomials of degree L can be obtained as a linear combination of the monomial 
functions (|161|) . and n\ + + ^3 = L. 

We also take into account constraint (|132|) via the following transformation: 

S (£1 + £2 + £3) = / 7^2 ex P W ■ (*i + £2 + £3)} • (162) 
2.7.2 Overlap matrix elements 

To calculate the overlap matrix B we consider the scalar product of two monomial 
functions ^{n},]* an d ^{ n '},k : 

<W I {"'}> = j dK f dp f f ^1^3 

X *{n},k*{n'},k exp {if • (£x + £ 2 + £3)} . (163) 
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Integration over R and Gaussian integration over p give a factor of 2ttA, where A 
is the area of the sample, and this yields 

3 

({n}\{n'}) = C 2 (27rA) ! J?L f[ M n . n ,{f), (164) 

J \ > j=l 



where 



M w (f) = J dt^'e-?'™** (165) 



We now express the complex variable £ in polar form as £ = £e l< ^ , and let 4> be the 
angle between f and £, i.e. 

<f> = fa - f . (166) 

Eq. (|165() then becomes 

/*27T POO 

M w (f) = e i(m - m,) ^ / #/ dec 1+mW 

Jo Jo 

x exp {i(m — m')4> + if£ cos ^ — £ 2 /6} . (167) 
The integration over c/> can be carried out using Bessel's integral |4()j : 



2tt 







e ±i(m-m')0+i/5cos0 = 27™1 W - m 'lj| m _ m ,[(/£), (168) 



where J m {z) is a Bessel function of order m. This simplifies Eq. 1)167(1 to 

M w (f) = 2^l"'le l ( m - m ')^^ mm ,(/), (169) 

where we define 

/>oo 

M mm >(f)= ^C 1+m+m 'e-^ /6 J lm - m/ \(m- (170) 



The integral in Eq. (|17()j) may be expressed |33j in terms of a confluent hypergeo- 
metric function <&(/?, 7; z) as 

2(^)/2 3 ( m W)/2 + l r(max{m;m / } + 1) 
Mmra' / = i TTj * 

|m — m'|! 

x <3? (max{m, m'} + 1, \m — m'\ + 1; — t) , (171) 

where t = 3/ 2 /2 and max{m, m'} is the largest of the integers m and m! . Applying 
the Kummer transformation 

*09,7;-t)=e-**(7-&7;t), (172) 
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simplifies Eq. (|171|) further. The fact that (7 — (3) is always a non-positive integer: 

7 — (3 = \m — m'\ — max{m, m'}, (173) 

means that $(7 — (3, 7;/;) reduces to a polynomial, and <£(/?, 7;— t) is then just a 
polynomial in t multiplied by e~ l . 

Integrating over cfif yields the final form of the overlap matrix elements 

poo 3 

{{n}\{n'}) = C 2 A(2^f5 LU / <lff\{i"- "'>!,. . v i.n. (174) 

Jo j=1 

where L = ni + n2 + n3. We define the constant C so that for L = L' = the matrix 
element ({71} | {n'}) = (0|0) = 1, which gives 

C 2 = , ~ttt - (175) 
3(2tt) 3 A v ; 

It should be noted that matrix elements in terms of symmetric polynomials can 
be constructed as a linear combination of the monomial matrix elements (j!74j) . 

2.7.3 Anyon-anyon interaction 

The interaction between anyons is of the form 

^k^ + ^ + ^}' (176> 

Only matrix elements for the first term V\2 will be calculated here as the others 
follow by analogy. Taking the Fourier transform of Eq. (1 176(1 yields 

Vi2(Z 1 ,£ 2 ) = l I ^ ex P^-^i-^)}, (177) 
which gives a matrix element 

<M \v 12 \ {«'}> = ^/ dR /^/(^i/^ / #1^8 

x *{„},k*{n'},k exp {if ■ (£1 + £ 2 + £3)} exp {iq ■ - £ 2 )} . (178) 
Following the method of Section f2.7.2l gives 



({n}\V 12 \ {n'}> = n^-^l /df /dq 

i=i ^ 7 



x exp {70 + (ni - 77^) + 70_(n 2 - re 2 ) + #f ( n 3 - n 3 )} 

x Mu„i(/+)M*ni (f-)M n3n > (/), (179) 
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where M. mm > are as defined in Eq. I|17l)|) . f± = f ± q, and <p± are the phases of f±. 

We now eliminate <j)f by the change of variables (f) = (f>{ — <^> q and tp± = 4>± — 4>f ■ 
Substituting for 0f , (f> + and (ft- in Eq. (|179|) and integrating over <^> q , we obtain 

({n} \V l2 \ { n >}) = C2A ^ )2 6 LL , n^-^l Td// 

/•27T 

x / d<p exp |z0(Z/ — L')} exp — n' x ) + i^-(n 2 — n' 2 )} 

Jo 

x M ni< (f+)M n2n , 2 (f-)M n3nl3 (f), (180) 

where ^-t and f± can be expressed in terms of the variables of integration, /, q and 

(f), as 

= (/ ± qe-*) /f±, f 2 ± = f 2 + q 2 ± 2fq cos <t>. (181) 

It is evident from Eq. ()18U|) that the anyon-anyon interaction matrix has the same 
block diagonal structure as the overlap matrix. The integrand in Eq. (|180|) can be 
further reduced to a product of a polynomial in q, f and e ±l<?i , and an exponential 
factor exp(— 9f 2 /2 — 3q 2 ). This integral can therefore be evaluated analytically for 
any {n} and {n'}. For the case of L = L' = 0, using Eq. (|175j) for C 2 , the anyon- 
anyon interaction matrix element reduces to 

(0\Vaa\0) = (182) 

2.7.4 Anyon-hole interaction 

The anyon-hole interaction is of the form 

V ah = -o \] — r + I — r + I 1 . 183 

3 [ \r lh \ \r 2 h\ \r3h\ J 

where we define 

r jh = £ j - p + hz. (184) 
Considering only V\h, we take the Fourier transform: 

Vlh (p, £i) = ~\j V M exp {iq • ($i - p)} , (185) 

where 

V^g) = —e~ qh . (186) 
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The matrix elements are then 

({n} \V lh \ {n'}> = - l -jdKjdpj J ^ 1 j j^Vah(q) J ^ 2 d£ 3 

x *{n},k*{n'},k exp {if ■ fa + £ 2 + £ 3 )} exp {iq • fa - p)} . (187) 

Integration over R is straightforward and gives a factor of A. However, the integra- 
tion over p is more complicated: 

/ dp exp {-(p-d) 2 /2-iq-p} = 27rexp(-g 2 /2-id-q). (188) 

Then, following the procedure in Section [2.7.21 once again, we obtain 

<{n} \V lh \ {n'}> = y df y dq VU^e^/^d-q (189) 

x exp - + - L') - (n\ - n^)]} 

X A4n 1 ni(/+)-Mn 2 ni(/)-Mn 3 ^(/), 

where M. m m' are as defined in Eq. 1)1701) . f + = f + q, and cft + is the phase of f + . As 
before, L = m + ni + 723. 

We now eliminate c/>f by substituting (ft = (ftf — (ftq and ip = (ft+ — (ftf. Choosing 
the x-axis along d, the integration over </> q is then 

/>2tt 

/ d0 q exp {z0 q (L — L') — idgcos0 q } = 27r(— i)' L ~ L ' J\L-L'\(dq), (190) 
which gives for the matrix elements: 



<w \v lh \ {„'» = jj^-nji e ^/*-«*j (dg) 

i=i 70 



/■oo / i 27r 

x / df f deft exp — L') + i%ft{n\ — n^)} 

where the new variables V an d /+ can be expressed in terms of the variables of 
integration as 

e*+ = (f + qe-**) //+, / 2 = f + g 2 + 2/gcos^. (192) 



For d = 0, the Bessel function J\L-L'\(dq) entering the integrand of Eq. (J191|) is 
non-zero only if L = L' . This indicates that for the case of zero in-plane momentum 
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(k = d = 0), the anyon-hole interaction matrix has the same block diagonal structure 
as the anyon-anyon and overlap matrices. Therefore, states with different angular 
momentum L z decouple. 

Eq. (|191|) may be simplified further using the explicit form of A4 mm /(/) in 
Eq. (|171l) and subsequent results. This yields the following expression for the 
anyon-hole interaction matrix elements for basis functions in terms of symmetric 
polynomials: 

/•oo 

{L,M\V ah \L',M') = - / exp(-3g 2 /2 - qh)J\ L - LI \{dq)Q LM ,UM'{q) dq, (193) 

Jo 

where Qlm,l'M'(q) is a polynomial in q (the lowest order polynomial is (Joo,oo = 1). 
This expression was previously derived in |47| . Note that for zero 2DEG-hole sep- 
aration (h = 0), Eq. (|193j) further reduces to an expression in terms of elementary 
functions. For the case of k = d = and L = V = 0, the anyon-hole interaction 
matrix element is then 



(oiv^io) = - 




(194) 



2.8 Six-particle anyon exciton 
2.8.1 Formulation 

We now introduce an anyon exciton consisting of a hole and five anyons with charge 
— e/5. We also make a boson approximation so that the statistical factor a = 0. 
Our justification for this step is as follows. It was shown in [17] that for large values 
of h (which is required for the AEM to be valid) the statistical factor a becomes 
unimportant, and the results for a = were very similar to those for a = ±1/3. We 
would therefore also expect this to be true for a = ±l/iV, as this is even closer to 
zero. From now on we consider a boson approximation (a = 0) for anyon statistics. 
The explicit form of the exciton basis functions is 

*L,M,k (R, P, {&}) = exp {ik • R + iz ■ [R x p] /2 - (p - d) 2 /4} 

5 

x Pl,m(Ci, ,C 5 ) II exp{-|C i | 2 /20} , (195) 
i=i 

where M enumerates all linearly-independent symmetric polynomials of degree L. 

To construct the symmetric polynomials Pl.m we need to consider only the 
elementary symmetric polynomials 02, 03, 04 and 05. Note that a\ = because 
of constraint 1)132(1 . From Eq. (|156j) . we find that the number of possible ways 
of constructing a polynomial of degree L is therefore the coefficient of x L in the 
expansion 
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Table 2: Possible ways of constructing a symmetric polynomial Pl from the elemen- 
tary symmetric polynomials <r 2 , 03, 04 and u 5 . 



Order, L 


No. of polynomials 


Structure 





1 


1 


1 







2 


1 


0"2 


3 


1 




4 


2 


n 

<T 2 , £7 4 


5 


2 


0"20"3> CT5 


6 


3 


of, cr 2 cr 4 , 0"| 


7 


3 


2 

0^03, CT2<75, 0"30"4 


8 


5 


4 2 2 2 
cj 2 , <T 2 (T4, cr 2 cr 3 , 0-305, cr 4 


9 


5 


of ^, 0"2°"5, CT2CT3CT4, CT3, <7 4 <7 5 


10 


7 


0"|, cr|a 4 , ofof , 0^030-5, cr 2 cr|, £7^(74, cr| 


11 


7 


O"!^, crl^S, 02C3O"4> <720"3, CT 2 (T4(T5, 
O3O5, 03O4 


12 


10 


fi 4 3 2 2 2 2 2 
O^, cr|o- 4 , cr^(7|, O^Os, crfcrf , O^CTgOM, 

O-20-f, Crf, CT 3 (T4<T5, O4 



• • t— 7 • — ^ = 1 + x 2 + x 3 + 2x 4 + 2x 5 

1 — x z 1 — x- 3 1 — x 4 1 — x° 

+ 3x 6 + 3x 7 + 5x 8 + 5x 9 + 7x 10 + ... . (196) 

The possible ways of constructing the first twelve polynomials are shown explic- 
itly in Table El It is evident that the number and complexity of the polynomials 
increases rapidly with the degree L. 

Following the method in Section 12.5.11 once again, we expand a general exciton 
wavefunction for given k in terms of the complete set of basis functions (|195|) : 

* = J><*i. (197) 

i 

As mentioned in Section 12.7.11 basis functions with different values of k and L are 
orthogonal, although basis functions with the same value of L but different M are 
not necessarily orthogonal, so we introduce the overlap matrix B of scalar products 
into the Schrodinger equation as follows: 



(198) 
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Note that B is block diagonal, and in the case of a six-particle exciton (see Table |2Jl, 
the block corresponding to L = 12 for example will be of size 10 x 10. As the 
Hamiltonian matrix H and the overlap matrix B are diagonal in k, we only need 
to consider the matrix elements in Eq. (|198|) that are diagonal in k. We use the 
following functions in monomials in the matrix element calculations, from which 
the basis functions in terms of symmetric polynomials can be obtained as a linear 
combination: 

*{n},k (R, P, {Q}) = Cexp {ik • R + iz ■ [R x p] 12- (p - d) 2 /4} 

5 

x CrCCCC II exp {-|0| 2 /20} , (199) 
i=i 

where {n} denotes the set of quantum numbers n\ to n§, for which it should be noted 
that n\ + ri2 + n% + 714 + n§ = L, and the constant C will be defined in Section r2.8.2l 
Constraint (|132ft is taken into account via the following transformation: 



*(*l+*2 + *3 + £4 + e5) = / 7^1 exp{if-(£i+£2 + £ 3 + £4 + £ 5 )}- (200) 



2.8.2 Overlap matrix elements 

To calculate the overlap matrix B, let us consider the scalar product of two monomial 
functions: 

<W I {"'}> = j dK j d pj / ^W«3^4^5 

X *{„},ktt{n'},k exp {if • + £ 2 + £ 3 + £4 + €5)} • ( 201 ) 

The integration over R and p gives a factor of 2ttA, and this leaves the following: 
(M I {*'}> = C 2 (2vrA) | 11^, (202) 



where 



(f) = y d£ ^ m C m 'e- 52/10+if ^. (203) 



If we now make the substitution £ = £e*^*, and let be the angle between f and 



i.e. 



0€ " </t, (204) 
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where fa is the phase of f , we obtain 

f2lT POO 



pzn poo 

M mm , (f) = exp {i(m - m')fa } / d<j> d£ £}+ m +™' 

Jo Jo 

x exp {i(m -m')</> + if £cos(/> - t? /10} . (205) 
We can now make use of Bessel's integral 

/ #e ±i(m - m,) * +i/5cos ^ = 2vril m - m 'lj| m „ m ,|(/0, (206) 
Jo 

to reduce our expression to the form 

M w (f) = 2^ m - m \^ m ~ m '^M mml (f), (207) 

where 

poo 

M mml {f)= dU 1+m+m 'e- e/W J lm ^\(fO- (208) 
Jo 

This can also be expressed in terms of a confluent hyper geometric function $(/?, 7; 2) 



as 



2 (^')/2 5 (>nW)/2 + l r(max{m)m / } + 1) _^ 
•M mTO '(j) = j 77j *' " 

|m — my. 

x $ (max{m, m'} + 1, \m — m'\ + 1; — t) , (209) 

where i = 5/ 2 /2 and max{m, m'} is the largest of the integers m and ml '. Eq. (|209j) 
can be simplified further if we apply the Kummer transformation: 

$(/3, 7 ;-t)=e-*$(7-A7;i) ) (210) 

and note that (7 — (3) is always a negative integer: 

7 — /3 = \m — m\ — max{m, m }. (211) 

This means that $(7 — /?, 7; t) reduces to a polynomial, and <&(/?, 7; — i) is then just 
a polynomial multiplied by e . 

If we now perform the integration over fa we obtain the final form of the overlap 
matrix elements 

poo 5 

({n}\{n>}} = C 2 A(2n) 5 5 LL , / df f \{ i'^M^' .(/), (212) 

- 70 7=1 
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where L = n% + n 2 + -\-ri4 + n^. Note that nothing depends on k in this expression 
for the matrix elements. It is convenient to define the constant C so that for L = 
L' = the matrix element ({n} \ {n'}) = (0|0) = 1. This yields 

° 2 - mWA- (213) 

In the above formulation we have only considered matrix elements in terms of 
monomial functions. As mentioned earlier, matrix elements in terms of symmetric 
polynomials can be easily constructed as a linear combination of the monomial 
matrix elements (j212j) . 

2.8.3 Anyon-anyon interaction 

The anyon-anyon interaction has the form 

^—sgi^ai- (214> 

We shall only calculate the matrix elements for the first term V12 as the others follow 
by analogy. We begin by taking the Fourier transform of Eq. 1)214(1 : 

Vi2(€i,&) = ^ / ^exp{iq -fo-fc)}. (215) 
This gives a matrix element in terms of monomials 

(WfelM)4/*/*/^/^ 

x exp {if • (^ + £ 2 + £ 3 + £ 4 + £ 5 )} exp {*q • - £2)} • (216) 
Following the procedure in Section [2.8.21 yields 

<W \V 12 \ {n>}) = n^-^l [dt /eiq 

x exp {i(j) + (ni — n[) + i0_(n 2 — n' 2 )} 

x exp {i0f [ (n 3 - ri 3 ) + (n 4 - n 4 ) + (n 5 - n' 5 ) ]} 

5 

x M nin/i (f + )M mnl2 (f-)llM n . n ,(f), (217) 

j=3 



X 
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where M. mm > are as defined in Eq. I|2U8|) . f-t = f ± q, and (f>± are the phases of f±. 

We now seek to eliminate <f>{ by the change of variables (j) = 4>f — c/> q and tp± = 
0± ~~ <t>f ■ After substitution and integration over ^> q we obtain 

({n}\V 12 \{n>}) = ^^5 LL ,fl^-^ Pdff Hdq P d* 
Zi) ~~[ Jo Jo Jo 

x exp {i<j)(L — L')} exp {iip + (ni — n[) + m/>_(ti2 — n'2)} 
x M ni „ ; (/ + ^„ 2 ^(/-)^^(/)K 4 „;(/)Ai„ 5 ^(/), (218) 
where ^ and f± can be expressed in terms of the variables of integration as 

e #± = ^ ± qe -i^ //±> /| = + g 2 ± 2/{? CQS ^ (219) 

Eq. (|218[) shows that the anyon-anyon interaction matrix has the same block diagonal 
structure as the overlap matrix. For the case of L = L' = 0, using Eq. (|213|) for C 2 , 
the anyon-anyon interaction matrix element reduces to 

(0\Vaa\0) = (220) 

2.8.4 Anyon-hole interaction 

The anyon hole interaction takes the form 

1(1 1 1 1 11 , . 

V ah = --{-. r + 1 r + : r + 1 f + 1 1 , (221) 

5 [ \r lh \ \r 2 h\ |r 3 /»| \r 4h \ \r 5h \ J 

where 

tjh = £j~P+ hz. (222) 
Considering only V±h, we take the Fourier transform: 

Vi h {p,£ 1 ) = -\ j j0^V ah (q)ex V {ic l -fa-p)}, (223) 

where 

V ah (q) = —e^ h . (224) 

q 

The matrix elements are then 



<{n} \V lh \ {n>}) = - 1 -JdRjdpJ J ^J j^VM 



x 



d£id€ 2 d£ 3 d€ 4 d£5 #{ n },k*{n'},k 



x exp {if ■ fa + £ 2 + £ 3 + £ 4 + | 5 )} exp {iq • fa - p)} . (225) 
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Integration over R and p gives 27rAexp (— q 2 /2 — id • q). Then, following the pro- 
cedure in Section 12,8.21 once again, we obtain 

({n}\V lh \{n'}) = - C2{2 f A f[^-^ fdff dqV 0/l (g)e-« 2 /^- q 

j= i J J 

x exp [i(j) + (ni - n[) + ifa [(L - L') - (rti - n[)] } 

5 

xM nini (/+)JJ^n^(/), (226) 

3=2 

where 7W mm ' are as defined in Eq. (|208jl . f + = f + q, and cp + is the phase of f + . As 
before, L = m + n2 + + 77,4 + 715. 

We now eliminate <f>f by substituting eft = <f>f — </> q and ip = 4> + — 4>f. Choosing 
d along the x-axis, the integration over </> q is then 

/ # q exp {i0 q (L - L') - idqeos^} = 2tt(-i)\ l ~ l 1 J\ L -L>\{dq), (227) 
•/ o 

which gives for the matrix elements: 

(m i^i {„'» = _^)^ H )i^'i n i'^i 

5 3=1 

f'OO f'QO 

x / d// / dqe^^ h J {L -L'\(dq) 
Jo Jo 

[■2lT 

x / d(f> exp {icp(L — L') + iip(ni — n^)} 
Jo 

5 

x^„ 1 <(/+)n- M n,^(/)' (228) 

3=2 

where -0 and /+ can be expressed in terms of the variables of integration as 

e**= (f + f 2 + = f 2 + q 2 + 2fq C o S <f). (229) 

Eq. (|228|) may be simplified further, and this yields the following expression for 
the anyon-hole interaction matrix elements for basis functions in terms of symmetric 
polynomials: 

POO 

{L,M\V ah \L',M') = - / exp(-5</ 2 /2 - qh)J\ L _ v Adq)Q LM ,L>M>(q) dq, (230) 
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where Qlm,l'M'(q) is a polynomial in q (the lowest order polynomial is (Joo,oo = 
1). For h = 0, Eq. ()23U|) further reduces to an expression in terms of elementary 
functions. For the case of L = V = 0, the anyon-hole interaction matrix element is 
then 



(o\v ah \o) = - 




(231) 



2.9 Exact results for k = 0, L = 

All the above results are simplified significantly for the state with k = and L = 0. 
For such a state, the (N + l)-particle wavefunction (|154j) with a = reduces to the 
form 

* (R, p, {d}) = C N exp {iz • [R x p] /2 - p 2 /4} exp j - £ |C P | 2 /4iV }> , (232) 



where the normalisation constant CV can be determined from the integral 
dfd^f^j dpe^^J d^-.-dtt 



IN 



exp <( £ (<€ p • f - £/2iV) } = ^(2^iV) W = 1. (233) 
p=i 



For the case of k = and L = 0, the interaction matrix elements of a (iV + 1)- 
particle exciton are straightforward to evaluate. The anyon-anyon matrix element 
is 

x J j^y ex P jE(*p-f-$/2W)+*q-(£i-& 

= i^l./Z. ( 234 ) 
4iV ViV V ; 



For non-zero inter-plane separation /i, the anyon-hole matrix element can be reduced 
to 
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A(2irN) 



X / (2^ 6XP |S (^P • f " ^/ 2iV ) + ^ ■ & -P)~* h j 

= ~[dq e -»*l*-* = "/^^erfc (h/V2N) , (235) 

where erfc(x) is the complementary error function. (Note that for h = 0, Eq. ()235|) 
can be simplified further to a numerical value.) Using the asymptotic expansion of 
erfc(x) it can be easily seen from Eq. ()235|) that (0|V^ft|0) — ► — as /i — ► oo, as 
expected. 

Eqs. (|234|) and (|235() also allow us to calculate the critical inter-plane separation 
h r at which the k = 0, L = state becomes unbound, i.e. when 



7T T 2 „ . . (JV - 1)V2 AvT 

- erfc(* c ) = [ -^J^, (230) 



where x c = h c /V2N. So, for iV = 3 the critical separation h c 5.39 I, for N = 5 
we find that /i c ~ 5.59 /, and for iV > 1 we have h c ~ 1.32\/2iV /. Notably, these 
critical separations are well inside the region for which the AEM is applicable. It 
should be emphasised that the state with L = is not the ground state for the 
anyon exciton at large separation h. For example, for a four-particle exciton 47 , 
the ground states for large separation satisfy a superselection rule L = 3m, where m 
is an integer, and when h — ► oo the ground state energy tends to its classical value, 
e cl = -(2/3f/ 2 /h. 

The model is applicable to real physical situations at large layer separations 
(h > I), although it remains solvable for all values of h, including h = 0. Moreover, 
it has been shown in [17] that the ground state of the four-particle problem (N = 3) 
at h = is the state with k = and L = 0. We expect the same to be true 
for N ^ 5, since the anyon-hole attraction will always overcome the anyon- anyon 
repulsion at small inter-particle separations. It can be shown (in a similar way to 
that in \47\ ) that the smallest average inter-particle separation corresponds to the 
L = case. Therefore, the case of k = 0, L = is of special interest since it predicts 
the anyon exciton binding energies in the limit h — ► 0. 

We are now in a position to write down a general expression for the binding 
energy of a (JV + l)-particle exciton at zero 2DEG-hole separation (h = 0): 



E b = - (OlKfclO) + (OIFJO) 



1 (N-l) 
2V2N 



w (237> 
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Eq. 1)237(1 has been written in this particular form to emphasise the key result that 
for any value of N there always exists at least one bound state of a neutral (N + 
l)-particle anyon exciton. For N = 1, Eq. (|237() yields the value which 
corresponds to that obtained in |1U1| for a standard diamagnetic exciton. For a 
four-particle exciton (N = 3), Ej, = (l — \/2/6) a/vt/6, which agrees with the result 
of [IT]. Finally, for N = 5 we have a binding energy of (1 - y/2/5)y/ir/10, which 
can also be obtained from Eqs. (|22()|) and (|231|) . 



2.10 Conclusions 

The anyon exciton model has been generalised to the case of a neutral exciton 
consisting of a valence hole and an arbitrary number N of anyons, and several 
important mathematical results have been obtained. A complete set of exciton basis 
functions has been obtained and these functions have been fully classified using a 
result from the theory of partitions. We have derived expressions for the overlap 
and interaction matrix elements of a six-particle system (N = 5), which describes 
an exciton against the background of an IQL with filling factor v = 2/5. In the 
particular case of k = and L = 0, we have found an expression for the binding 
energy of a (A r + l)-particle exciton at zero anyon-hole separation, which agrees with 
known results for a standard diamagnetic exciton and four-particle anyon exciton. 
We have also shown that the {N + l)-particle exciton remains bound for 2DEG- 
hole separations exceeding several magnetic lengths, when the anyon exciton model 
becomes applicable to real physical systems. 

3 Appendices 

3.1 Solution of real-space Schrodinger equation in two dimensions 

We apply the method of separation of variables to Eq. (|39|) , making the substitution 

(p) = R(pM<i>). (238) 
Introducing a separation constant m 2 , we can obtain the angular equation 



d 2 ^ 

dcj) 2 

with the solution 



+ m 2 <$> = 0, (239) 



= -^e im *. (240) 



The corresponding radial equation (with E = — qfi) is 



d 2 R ldR (2 2 m 2 , r , , , 
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We make the substitution 

R(p) = Cp\ ml e- qop w(p), (242) 
where C is a normalisation constant. This leads to the equation 

p^ + (2\m\ + 1 - 2q p)^ + (2 - 2\m\q - q )w = 0. (243) 
Making a final change of variables (5 = 2qop, we obtain 

/3 ^ + (2W + 1 _ /3) ^ + (I_ W _l)„ = 0. (244) 

This is the confluent hyper geometric equation [33], which has two linearly inde- 
pendent solutions. If we choose the solution which is regular at the origin, then 
this becomes a polynomial of finite degree if qo = (n + 1/2) -1 , with n = 0, 1, 2, . . . 
Eq. (|244|) then becomes the associated Laguerre equation 4Q] , the solutions of which 
are the associated Laguerre polynomials 

^ = ^'i|(/3)=^ l HL|(2M- (245) 
We can now write the real-space wavefunction in the form 

*W(p) = ^p^e-^L 2 ^L(2q p)e im ^, (246) 



where the reason for the subscript on (ft is explained in Section [1.7. 31 

To normalise this wavefunction we need to make use of the integral |4(J| : 



e- 2 «^(2gop) 2|m|+1 L^ | (2g p)L^L|( 2 9oP) d(2q oP ) 

(n + \m\)\ 
(n — \m\)\ 

The normalised wavef unctions are therefore: 



(2n + l). (247) 



= ftn + H)! (^P) |m| ^lt|(2 g 0P)e^, (248) 
satisfying the following orthogonality condition: 

J K 1 m 1 (p) q 'n 2m2 (p)dp = 5 nin2 5 mim2 . (249) 
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3.2 Derivation of expression for A 2 
From Eq. (|41|) we have 

~ - 2 I 2 
(q x L z - L z x q) p 

P . 



(250) 



= [2(q x L z ) — iq\ 2 — -p ■ [2(q X L z ) - iq] - - [2(q X t z ) - tq] • p + 4. 

We further expand as follows: 

[2(q x t z ) - iq] 2 = 4(q x t z f - 2iq • (q x t z ) - 2t(q x L z ) • q - q 2 (251) 
= 4q 2 L 2 + 2 q 2 - q 2 = q 2 (4L 2 + l), 

and 

--p • [2(q x t z ) - iq\ - -[2(q x L z ) - jq] • p = --(4L 2 + 1). (252) 
P P P 

Substituting Eqs. (|23T|) and into Eq. (j23Uj) gives 

which, from Eq. is just 

A 2 = F(4L 2 + 1) + 4. (254) 
3.3 Variable phase method 

In this Section, which is based on the formulation in we derive the basic equa- 
tions of the variable phase approach in two dimensions from the radial Schrodinger 
equation. This derivation is similar to that in three dimensions jll)5j . 

3.3.1 Scattering phase shifts 

The relative in-plane motion of two interacting particles with masses M a and 
and energy of relative motion E can be considered as the motion of a particle 
with mass \x a \, = M a Mb/(M a + M^) and energy E, moving in an external central 
potential V(p). This motion is described by a wavefunction satisfying the stationary 
Schrodinger equation 

h 2 /l d d 1 d 2 \ 
2/iafe \pdp dp p 2 dp 2 J 
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Owing to the axial symmetry of the potential V(p), we can separate variables in the 
expression for the wavefunction 

if> m (p,<p) =R m (p)e imip , m = 0,±l,±2,... . (256) 
The equation for the radial function R m (p) reads 

R n m + - p R' m + (k 2 - U(p) -^R m = 0, (257) 

where k 2 = 2p a i ) E/h 2 and U(p) = 2/j, a i,V(p)/h 2 . In what follows we consider m > 
only, as R- m (p) = R m (p). 

We assume that the interaction potential vanishes at infinity (the precise decay 
rate will be discussed later). Then, at large distances the radial function satisfies 
the free Bessel equation, whose general solution is 

R m (p) = A m [J m (kp) cos5 m - N m (kp) sm5 m ] (258) 

oo / 2 \ 1/2 

A m — cos(kp- (2m + l)vr/4 + (5 m ), 

\TTKp J 

where S m is the scattering phase shift and J m (kp) and N m (kp) are the Bessel 

and Neumann functions, respectively. 

In the variable phase approach, A m and 8 m are considered not as constants, 
but as functions of the distance p. The amplitude function A m {p) and the phase 
function 5 m (p) are introduced by the equation 

Rm(p) = A m {p)[J m (kp) cos 5 m (p) - N m (kp) sin 5 m (p)\, (259) 

with an additional condition, which we are free to choose as 

R' m (p) = A m {p)[j' m (kp) cos 5 m (p) - N' m (kp) sin S m (p)], (260) 

where the prime indicates differentiation with respect to p. The phase function 5 m (p) 
has a natural physical interpretation as the phase shift produced by a potential cut- 
off at a distance p. 

Differentiating Eq. (|26()|) and substituting the resulting expression, together with 
Eqs. M and M . into Eq. (12371) . we get 

A' m (p)[J' m {kp) cos 5 m (p) - N' m {kp) sin 8 m (p)] (261) 
- 5' m {p)A m (p)[J' m (kp) sin 5 m {p) + N' m {kp) cos 5 m (p)} 
= U(p)A m (p)[J m (kp) cos5 m (p) - N m (kp) smS m (p)]. 

To obtain Eq. (|261|) we used the fact that the functions J m (kp) and N m (kp) satisfy 
the free Bessel equation 

K + -F' m + (V - F m = 0. (262) 
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Equating the derivative of Eq. (|259j) to Eq. (|260j) implies the following condition 
on the derivatives of the amplitude and the phase functions: 

A 'm(p)[ J m(kp) COS 6 m (p) ~ N m (kp) SUl8 m (p)] 

= 8' m (p)A m (p)[J m {kp) sin 8 m (p) + N m (kp) cos 8 m (p)]. (263) 
Substituting A' m (p), obtained from Eq. ()263|l . into Eq. 1)261 Jl yields 

- 8' m {p)[J m {kp)N' m {kp) - N m (kp)j' m {kp)} 

= U(p)[J m (kp) cos 8 m (p) - N m (kp) sin 8 m (p)} 2 . (264) 

Eq. (|264|) can be simplified further, using the Wronskian of the Bessel functions 

W{ J m (x), N m (x)} = J m (x)^-N m (x) - N m (x)^-J m (x) = — , (265) 

dX CIX 7TX 

and thus becomes 

-r8 m {p) = ~ pU(p)[J m (kp) cos 8 m (p) - N m (kp) sin 5 m (p)} 2 . (266) 
dp 2 

This phase equation, Eq. (|266|) . is a first-order nonlinear differential equation of the 
Ricatti type, which must be solved with the initial condition 

MO) = 0, (267) 

thus ensuring that the radial function does not diverge at p = 0. The total scattering 
phase shift 8 m can be obtained as a large-distance limit of the phase function 5 m (p): 

8 m = lim 8 m (p). (268) 

For numerical convenience, instead of the initial condition Eq. I|267j) . the small-/? 
expansion is used: 

MP) - ~ 22m+1(m!)2 I U(p') P ' 2m+1 dp\ P - 0. (269) 

From Eq. ()266|) and the asymptotic expansions of the Bessel functions, one can 
see that the variable phase method is applicable only if the scattering potential U (p) 
satisfies the necessary conditions 

POO 

/ U{p')dp' -»■ 0, p -» oo, (270) 
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and 

f?U(p)^0, p^O. (271) 

The Stern-Howard potential V s (p), defined by Eq. (^Q), behaves like at small 
distances and like p~ 3 at large distances. Such behaviour allows the application of 
the variable phase method to this potential. 

3.3.2 Bound-state energies 

For states with negative energy of the relative motion (bound states), the wavenum- 
ber k is imaginary {k = in), and we introduce the function r] m (p, k) which vanishes 
at the origin and satisfies the nonlinear equation 



j r r, 1 2 

d , s vr 



— 7]m(p,K) = --pU(p) 



2 

I m (np) cos rfm(p, «) + -K m (Kp) sin rj m (p, k) 

7T 



(272) 



where I m {np) and K m {np) are modified Bessel functions of the first and second kind, 
respectively. Eq. Q272|) is derived in the same fashion as Eq. 1)266(1 . The functions 
I m (np) and K m {K,p) represent two linearly-independent solutions of the free radial- 
wave Schrodinger equation for the negative value of energy E = —h 2 n 2 /2p a i ) , and 
cot 77 m characterises the weights of the diverging [I m (Kp)] and converging [K m (np)] 
solutions as p — > oo. For the bound state the diverging solution vanishes, implying 
the asymptotic condition 

r] m (p ->■ oo, k v ) = (v - 7T, v = 1,2,3, ... . (273) 

Here, v numerates the bound states for a given m and (y — 1) is the number of non- 
zero nodes of the radial wavefunction. For numerical solution of Eq. 1)272)1 . instead 
of the boundary condition r) m (Q,n) = 0, an asymptotic initial condition (analogous 
to the condition Eq. ()269)1 for the phase function 5 m (p)) is used. 
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